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Abstract 

In this paper we present an elementary derivation of the semi-classical spec- 
trum of neutral particles in a field theory with kink excitations. In the non- 
integrable cases, we show that each vacuum state cannot generically support 
more than two stable particles, since all other neutral exitations are reso- 
nances, which will eventually decay. A phase space estimate of these decay 
rates is also given. This shows that there may be a window of values of 
the coupling constant where a particle with higher mass is more stable than 
the one with lower mass. We also discuss the crossing symmetry properties 
of the semiclassical form factors and the possibility of extracting the elastic 
part of the kink S-matrix below their inelastic threshold. We present the 
analysis of theories with symmetric and asymmetric wells, as well as of those 
with symmetric or asymmetric kinks. Illustrative examples of such theories 
are provided, among others, by the Tricritical Ising Ising, the Double Sine 
Gordon model and by a class of potentials recently introduced by Bazeira et 
al. 



On leave of absence from International School for Advanced Studies, Trieste (Italy) 



1 Introduction 



The aim of this paper is to explain in the easiest possible way the presence of neutral 
bound states in two dimensional field theories with kink topological excitations. The 
theories that we will consider are those described by a scalar real field (f{x), with a 
Lagrangian density 

^ = lid,^f-U{^) , (1.1) 

where the potential U{ip) possesses several degenerate minima at ip^a^ (a = 1, 2, . . . , n) 
as the one shown in Figure 1. These minima correspond to the different vacua | a) 
of the associate quantum field theory. 




Figure 1: Potential U{(f) of a quantum field theory with kink excitations. 

The basic excitations of this kind of models are kinks and anti-kinks, i.e. topological 
configurations which interpolate between two neighbouring vacua. Semiclassically 
they correspond to the static solutions of the equation of motion, i.e. 

dl^{x) = Uy{x)] , (1.2) 

with boundary conditions (p{—oo) = ip^a^ and ip{+oo) = where 6 = a ± 1. This 
equation can be equivalently expressed in terms of a first order differential equation 

dip 



dx ±v/2Z7M . (1.3) 

Denoting by (pab{x) the solutions of this equation, their classical energy density is 
given by 

^ab{x) = I (^] + U{ipab{x)) , (1.4) 



2 \ dx 
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and its integral provides the classical expression of the kink masses 

/oo 
eab{x) . (1.5) 
-oo 

As a rule of thumb, it is useful to notice that the classical masses of the kinks (pab{x) 
are simply proportional to the heights of the potential between the two minima (f'^a^ 
and ipf\ 

The classical solutions can be set in motion by a Lorentz transformation, i.e. 
fab{x) '^ab [{x ± vt) /Vl — f^J • In the quantum theory, these configurations de- 
scribe the kink states | Kab{d) ), where a and h are the indices of the initial and final 
vacuum, respectively. The quantity 9 is the rapidity variable which parameterises 
the relativistic dispersion relation of these excitations, i.e. 

E = MabCos\ie , P = M„6sinh^. (1.6) 

Conventionally | Ka,a+i{,G) ) denotes the fcmfc between the pair of vacua {| a ), | a + 1 )} 
while I Ka+i,a ) is the corresponding anti-kink. For the kink configurations it may 
be useful to adopt the simplified graphical form shown in Figure 12^. 

I n > 

I a+l> I 

I a > 1 

I0> 

^ a,a+l ^ a+l,a 



Figure 2: Kink and antikink configurations. 

The multi-particle states are given by a string of these excitations, with the adja- 
cency condition of the consecutive indices for the continuity of the field configuration 

|i^a„a.(^l)i^a.,a3(^2 )i^a3.a4(^3)---) , (a.+l = ± 1) (1-7) 

^Although very convenient, one should keep in mind that this graphical representation over- 
simplifies the exponential approaching to the vacua, which can be different for the one to the right 
and for the one to the left. 
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In addition to the kinks, in the quantum theory there may exist other excitations in 
the guise of ordinary scalar particles (breathers). These are the neutral excitations 
I Bc{6) )a (c = 1,2, . . .) around each of the vacua | a). For a theory based on a 
Lagrangian of a single real field, these states are all non-degenerate: in fact, there 
are no extra quantities which commute with the Hamiltonian and that can give 
rise to a multiplicity of them. The only exact (alias, unbroken) symmetries for a 
Lagrangian as (jl.ip may be the discrete ones, like the parity transformation P, for 
instance, or the charge conjugation C. However, since they are neutral excitations, 
they will be either even or odd eigenvectors of C. 

The neutral particles must be identified as the bound states of the kink-antikink 
configurations that start and end at the same vacuum | a), i.e. | Kab{Oi) Kba{G2) ), 
with the "tooth" shapes shown in Figure El 



I0> ' ' 1 I 

a 

I > I 

a-1 

Figure 3: Kink-antikink configurations which may give rise to a bound state nearby 
the vacuum | )„. 

If such kink states have a pole at an imaginary value i within the physical strip 
< Im^ < vr of their rapidity difference 6 = 6i — 62, then their bound states 
are defined through the factorization formula which holds in the vicinity of this 
singularity 

I Kab{9,)KU02)) ^ ^TT^^ \Bc)a . (1.8) 

In this expression (7^^ is the on-shell 3-particle coupling between the kinks and the 
neutral particle. Moreover, the mass of the bound states is simply obtained by sub- 
stituing the resonance value i within the expression of the Mandelstam variable 
s of the two-kink channel 

9 

s = 4M2, cosh^- m, = 2MabCos^. (1.9) 

Since the kink (pa-i,aix) that interpolates to the vacuum on the left of | a) may 
be different from the kink ipa,a+i{x) which interpolates to the vacuum on the right. 
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a-priori there could be two different towers of breathers that pile up in each vacuum: 
the first coming from the poles of | -ft'a,a-i(^i)-^a-i,a(^) ) while the second one from 
the poles of | Ka^a+i{d)Ka+i^a{d2)) ■ Howevcr, as we will discuss in the following, this 
situation cannot occur. 

Concerning the vacua themselves, as well known, in the infinite volume their clas- 
sical degeneracy is removed by selecting one of them, say \ k), out of the n available. 
This happens through the usual spontaneously symmetry breaking mechanism, even 
though - stricly speaking - there may be no internal symmetry to break at all. This 
is the case, for instance, of the potential shown in Figure 1, which does not have 
any particular invariance. In the absence of a symmetry which connects the various 
vacua, the world - as seen by each of them - may appear very different: they can 
have, indeed, different particle contents. The problem we would like to examine in 
this paper concerns the neutral excitations around each vacuum, in particular the 
question of the existence of such particles and of the value of their masses. 

The answer provided by the perturbation theory to this question is straightfor- 
ward: after defining t]{x) = — ^pf^ and making a Taylor expansion of U{ip) near 

U{^^^^+V) = \ujlr,^ + \\,r,' + \\Ari' + --- (1.10) 

one identifies the mass m of the fundamental particle around the vacuum | a ) and 
ujai while the rest of the expansion with its interaction terms. The quantity Ua 
is of course the zero-order value of m of such a particle, but the crucial point is 
another one: according to the perturbation theory, as far as the potential has a 
quadratic curvature at its minimum, there is always a neutral excitation above the 
corresponding vacuum state. This conclusion is, unfortunately, false. 

A famous counter-example is given by the Sine-Gordon model, i.e. the quantum 
field theory associated to the potential 

2 

Usg{^) = [1 - cos(/?(^)] , (1.11) 

where mo is a mass-like parameter and is a coupling constant. Such a theory 
has an infinite number of degenerate vacua \ a) {a = 0, ±1, ±2, . . .), localised at 
(f^^ = 27Ta/(3, each of them with the same curvature u"^ = rn^. Through the above 
perturbative argument, one would conclude that each vacuum has always, at least, 
one neutral excitation. On the other hand, the exact S'-matrix of this model P 
shows that the situation is rather different: indeed, such a particle does not exist if 
/32 > An. 

As in the case of the Sine-Gordon model, the knowledge of the exact S'-matrix of 
a quantum problem would obviously provide a clear cut answer to the question of the 
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particle content of a theory: a proper identification of its poles gives its spectrum. 
This has been amply proved by the large number of the exact S'-matrices associated 
to the integrable deformations of Conformal Field Theories [2] (for a review, see j3]). 
But, what happens if the theory is not integrable? How can we proceed if the exact 
S matrix is not known? 

2 A semiclassical formula 

The particle content of certain non-integrable models can be studied by using the 
so-called Form Factor Perturbation Theory IHEHEl- As shown in this approach 
can be also extended to compute, in a reliable way, the decay widths of the unstable 
particles of the theory^. In addition to this approach, another interesting route for 
investigating the non-integrable models comes from semi-classical methods. Origi- 
nally proposed by Dashen-Hasslacher-Neveu P and by Goldstone-Jackiw (TU], this 
approach has been recently applied either to study quantum field theories on a finite 
volume ^21 ^1 to obtain their spectrum at the semiclassical level ^3] . 

The main difference between the two approaches is the following. The Form 
Factor Perturbation Theory is a formalism based on the S'-matrix theory. More 
precisely, it moves its first steps with the exact scattering amplitudes of the inte- 
grable models, reached as a limit of the non-integrable ones. On the contrary, the 
semi-classical methods build their analysis on the Lagrangian density of the model, 
irrespectively whether it describes an integrable system or not. As shown in ^3], 
the two approaches sometimes coincide while, in some other cases, they complement 
each other, i.e. one needs both methods to recover the whole spectrum of the theory. 

For the problem that concerns this paper - to find the neutral spectrum of 
the Lagrangian theory (jl.ll) - the semiclassical methods are the obvious choice. 
Their correct implementation may need, though, important pieces of information 
coming from the S'-matrix theory or from the Form Factor Perturbation Theory. 
The starting point of our analysis is a remarkably simple formula due to Goldstone- 
Jackiw ^ni- In its refined version, given in jTTj and rediscovered in jT^j, it reads as 

^The analitic prediction of the paper [7] for the decay widths of the unstable particles of the 
Ising model has been confirmed by the numerical analysis done in [HI. 
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follows^ (Figure IH) 

/oo 
cia;e*^^-^>,b(x) , (2.12) 
-oo 

where 6 = 6i — 62. 




K 

ab 

Figure 4: Matrix element between kink states. 

Notice that, if we substitute in the above formula 6 iir — 9, the corresponding 
expression may be interpreted as the following Form Factor 

Kbi^) = /(^TT-^) = (a I ^(0) I KabiOi) KUO2)) . (2.13) 

In this matrix element, it appears the neutral kink states around the vacuum | a) 
we are interested in. 

By following the references [TUllTTllT^llSj . let's firstly recall the main steps that 
lead to this formula and let's make the first comments on its content. Denoting 
the adimensional coupling constant of the theory generically by g, we will assume 
that the mass of the kink becomes arbitrarly large when g ^ 0, say as Mat — 1/g. 
Consider now the Heisenberg equation of motion satisfied by the field (p{x) 

(d^-dl) ^ix,t) = -U'[cfix,t)] , (2.14) 

and sandwich it between the kink states of momentum pi and p2. By using 

{KM I if{x,t) I KM) = e-'^P'-P'^-^" {KM I m I KM) (2-15) 

■^The matrix element of the field ip{y) is easily obtained by using ^p{y) ~ e^^^^*^^ (y9(0) e'^''*''' 
and by acting with the conserved energy-momentum operator on the kink state. Moreover, for 
the semiclassical matrix element F^f^{9) of the operator ^^[(^(O)], one should employ G[ipabix)]. For 
instance, the matrix element of <^^(0) are given by the Fourier transform of (/9^j(a;). 
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we have 



[-iPl-P2)APl-P2r] {KabiPl) \m \Kabip2)) 
= -{KabiPl) I U'MO)] I Kab{p2)) . (2.16) 

Once it has been extracted the x^-dependence of the matrix element (j2.15p . the 
remaining expression {Kabipij \ V5(0) | Kab{p2) ) should depend on the relativistic 
invariants of the channel of the two kinks. Since these invariants can be expressed 
in terms of difference of the rapidities of the two kinks, this suggests to adopt the 
rapidity variables and write eq. ()2.1(i|l as 

2Mf,(cosh^-l) {Kab{9i) |</^(0) I Kab{e2)) 
= -{KaM I U'l^fiO)] I KaM) , (2.17) 

where 6 = 6i — 62- Let's now assume that the kinks are sufficiently slow, so that 
their dispersion relations can be approximated by the non-relativistic expressions 

E = M coshe ~ M 1 + — , P = M sinh^ ^ < M . (2.18) 



In this quasi-static regime, we can define the matrix element of the field (/?(0) be- 
tween the kink states as the Fourier transform with respect to the Lorentz invariant 
difference 6 = 61 — 62 

/oo 
dx e'''^^'^ fix) , (2.19) 
-00 

with the inverse Fourier transform given by 

m = I ^e-^^^-^^7r.W • (2.20) 
In the quasi-static limit, the left hand side of equation ()2.17j) becomes 

Mle'fm , (2.21) 
which, in real space, corresponds to 

'^'^ fabix) . (2.22) 



dx'^ 

Concerning the right hand side, let's assume that f/'[0(O)] can be expressed in terms 
of powers of the field (p{x) (either as a finite sum or an infinite series) 

U'[ip\ = ^a„<^" . (2.23) 

n=l 
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Consider now the matrix elements of the generic term of this expression between 
the kink states 

{Kab{e^) I <^"(0) I Kabie^)) . (2.24) 
By inserting {n — \) times a complete set of state, we have 

^ (/Cfe(^i) I ^(0) I mi) (mi I ^(0) I ma) . . . (m„_i | ^(0) | Kab^e^)) . 

mi,...m„_i 

(2.25) 

The only states which are involved in the above sums are those having the same 
topological charge of the kink Kab-, with the lowest mass states given precisely by 
the kinks Kab{Oi) themselves. By truncating the sums just on these states and using 
the definition (|2.19p . we have then 

rfxe^"^-^M/(x)) . (2.26) 

oo 

Hence, at the leading order 1/g, the function /(x) satisfies the same differential 
equation ()1.2|1 satisfied by the static kink solution, i.e. 

-^U^) = U'lfabix)] , (2.27) 

arriving then to the result ()2.12j) . 

The appealing aspect of the formula ()2.12|) stays in the relation between the 
Fourier transform of the classical configuration of the kink, - i.e. the solution ipab{x) 
of the differential equation p.Hj) - to the quantum matrix element of the field ^p{0) 
between the vacuum | a) and the 2-particle kink state | Kab{Oi) -ft'ba(^2) )• Once the 
solution of eq. ()1.3|) has been found and its Fourier transform has been taken, the 
poles of Fab{0) within the physical strip of 9 identify the neutral bound states which 
couple to if. Then, their mass can be extracted by using eq. (jl.9j) . while the on-shell 
3-particle coupling g^^^ can be obtained from the residue at these poles (FiguraE)) 

hm^ {e-iul,)F,,{e) = igl, {a \ ^(0) \ B,) . (2.28) 

ab 

It is important to stress that, for a generic theory, the classical kink configuration 
9?ab(a;) is not related in a simple way to the anti-kink configuration ipi,a{x). It is 
precisely for this reason that neighbouring vacua may have a different spectrum of 
neutral excitations, as shown in the examples discussed in the following sections. 

It is also worth noting that this procedure for extracting the bound states masses 
permits in many cases to avoid the semiclassical quantization of the breather solu- 
tions in], making their derivation much simpler. The reason is that, the classical 
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Figure 5: Residue equation for the matrix element on the kink states. 

breather configurations depend also on time and have, in general, a more complicated 
structure than the kink ones. Yet, in non-integrable theories these configurations 
do not exist as exact solutions of the partial differential equations of the field theory. 
On the contrary, in order to apply eq. ()2.12j) . one simply needs the solution of an 
ordinary differential equation, the one given by in the absence of an exact 

expression, one could even conceive the idea of employing the solution extracted 
by a numerical integration of the equation ()1.3|) . an operation which requires few 
seconds on any laptop"^. 

Let's now add a remark of technical nature: the Fourier transform (j3.32|) has 
always a singular part, due to the constant asymptotic behaviours of at x ^ 
±00. This piece can be easily isolated by splitting the integral as follows 



fab{k) = I dxe"'''(pab{x) = / dx e'""^ ifabix) + dx e'"'' ifabix) 

J -00 Jo 

POO 

dx e'"^ [ipabix) - <^i°)] + / dx e'"^ ipabix) - iff^ (2.29) 

Jo L ^ 

/O poo 
dxe'^'^' + ^f / dxe'^"" . 
■00 Jo 

The singular part is enconded in the last two terms, for which we have 

fO 







''^One might doubt about the possibihty to locate the poles of Fab{d) by performing a Fourier 
transform of the numerical solution. However, their position can be easily determined in a different 
way, i.e. by looking at the exponential behavior of the solutions at x — *■ ±00, as discussed in the next 
section. This behavior can be quickly extracted by numerical solutions and it is also, analitically, 
perfectly under control. 
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where P stays for the Cauchy principal value. Since these singular terms do not 
contain information on the pole structure of the matrix elements, they will always 
be discarded from now on. Concering the regular part, it can be computed by means 



Additional comments on the nature of the poles of the semiclassical form factors 
will be given in the next section, with the help of an explicit example. 

The range of validity of the formula ()2.12j) is a more delicate issue. As discussed 
above, its derivation relies on the basic hypothesis that the kink momentum is very 
small compared to its mass, and also on the possibility of neglecting intermediate 
higher particle contributions ^01 ^] • These two assumptions usually translate 
into the combined condition 6 ~ 0{g) ^ 1, where g is the adimensional coupling 
constant of the theory. This authorises, for instance, to substitute in the result of 
the Fourier transform, 6 sinh6' (since 6 is infinitesimal), but keeping untouched 
all expressions containing 6/g. But, the above constraint may result in a different 
level of accuracy on various physical quantities, with the precision that may also 
depend on the model under investigation. 

Consider, for instance, the Sine-Gordon model, a theory that will be studied 
in details in Section |S| Its semiclassical mass formula turns out to be an exact 
expression, i.e. valid at all orders in the coupling constant till its critical value 
f3^ = 4:71, beyond which, the bound states disappear. Moreover, the numerical values 
of the semiclassical Form Factors do not significatively differ from the exact ones for 
all real values of 6. On the wave of these success, one may dare to extract the S- 
matrix by using the semiclassical Form Factors. If one tries to do so, the outcoming 
expression turns out to be remarkably close to the exact one but, it always fails to 
meet an important crossing symmetric factor. The same happens in other models 
too, as it will be explained later. 

All this to say that, some caution is necessary in handling the results obtained 
by using the formula ()2.12j) . Morally it has the same status of the WKB approx- 
imation in quantum mechanics: this usually provides very accurate results for the 
discrete states, remaining nevertheless a poor approximation for those of the con- 
tinuum. But even for the bound states, the formula ()2.12|) seems sometimes to lead 
to some puzzling conclusions. This happens, for instance, in the case of the Double 
Sine-Gordon model fSlElEI- Fortunately, there is a way out from the possible 
paradox that involves the bound states, the solution of this problem being one of 
the main motivations of this paper. In summary, properly handled, the semiclassical 



of the derivative of the kink solutions, by using = —ikT{f), where 
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formalism remains one of the most powerful method for extracting the spectrum of 
the bound states in kink-like theories. 

In the next two sections we will first analyse a class of theories with only two 
vacua, which can be either symmetric or asymmetric ones. After the analysis of 
the Sine-Gordon model, we will proceed to discuss the interesting case of a vacuum 
state, in communication through asymmetric kinks, with two of its neighbouring 
ones. The simplest example of this kind of situation is provided by the Double Sine 
Gordon model. As we shall see, the conclusions drawn from all the above cases 
enable us to address the study of the most general theories with kink excitations, 
as shown by the examples discussed in the last section. 



3 Symmetric wells 



A prototype example of a potential with two symmetric wells is the theory in its 
broken phase. The potential is given in this case by 

4 V X 



U{^) = t(v'-'^] . (3.30) 



Let us denote with | ±1 ) the vacua corresponding to the classical minima f/^l^'* = 
±;^. By expanding around them, (p = (p^^^ + rj, we have 

f/((^i°^ + T]) = T]"^ ± mVXT]^ + ^r/^ . (3.31) 

Hence, perturbation theory predicts the existence of a neutral particle for each of 
the two vacua, with a bare mass given by = \/2m, irrespectively of the value of 
the coupling A. Let's see, instead, what is the result of the semiclassical analysis. 
The kink solutions are given in this case by^ 



mx 



a = ±1 (3.32) 



and their classical mass is 



2^/2 rrv^ 



Mo = / e(x) dx = ^ — . (3.33) 



oo 



The value of the potential at the origin, which gives the height of the barrier between 
the two vacua, can be expressed as 

f/(0) = l^^o , (3.34) 



^In the following we will always discard the integration constant xq on which is localised the 
kink solution. 
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and, as noticed in the introduction, is proportional to the classical mass of the kink. 

If we take into account the contribution of the small oscillations around the 
classical static configurations, the kink mass gets corrected as jH] 
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It is convenient to define 

c= ( — 1 >0 



2ti 4v^ 
and also the adimensional quantities 



3A , _ g 

2TTni? ' 1 — -Kcg 



9 = ; e = • (3.36) 



In terms of them, the mass of the kink can be expressed as 

M = = . (3.37) 

IT t, Tc q 

Since the kink and the anti-kink solutions are equal functions (up to a sign), their 
Fourier transforms have the same poles. Hence, the spectrum of the neutral particles 
will be the same on both vacua, in agreement with the Z2 symmetry of the model. 
For explicitly computing it, let's first consider 



2 



^1 = ^ , (3.38) 

dx)_^^^ ySAcosh^^ 

and then apply the prescription for the Fourier transform given in the previous 
section. As a result, we have 



By making now the analitical continuation 6 in— 6 and using the above definitions 
fl3.36|) . we arrive to 

F.aA^) = {a I ifiO) I K.a,a{Oi)Ka,-a{02)) = / x , (3.39) 



sinh 

where Aa is a constant, given by 

'3(l + vrcO^'/' 



(in-e) 



Aa = iam 



(^0 
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The poles of the above expression are located at 



Or, 



i-K {1 — S, n) 



n 



0,±1,±2, 



(3.40) 



and, if 

e > 1 , (3.41) 

none of them is in the physical strip < Im^ < tt. Consequently, in the range of 
the coupling constant 



A 2tx 1 

— > 

^2-31 



= 1.02338... (3.42) 

<j i + TTC 

the theory does not have any neutral bound states, neither on the vacuum to the 
right nor on the one to the left. Viceversa, if ^ < 1, there are n = j neutral 
bound states, where [x] denote the integer part of the number x. Their semiclassical 
masses are given by 



(n) 





■ 




sin 


n — 


= n nif, 




2 



\ 3 A2 2 

'-32;;^" + 



(3.43) 



Note that the leading term is given by multiples of the mass of the elementary boson 
\ Bi). Therefore the n-th breather may be considered as a loosely bound state of n of 
it, with the binding energy provided by the remaining terms of the above expansion. 
But, for the non-integrability of the theory, all particles with mass m„ > 2mi will 
eventually decay. It is easy to see that, if there are at most two particles in the 
spectrum, it is always valid the inequality 1112 < 2mi. However, if ,^ < |, for the 
higher particles one always has 



mfc > 2mi 



for = 3, 4, . . . n . 



(3.44) 



According to the semiclassical analysis, the spectrum of neutral particles of ip^ theory 
is then as follows: (i) if > 1, there are no neutral particles; (ii) if ^ < ^ < 1, there 
is one particle; (iii) if | < ^ < | there are two particles; (iv) if ^ < | there are | 
particles, although only the first two are stable, because the others are resonances. 



The decay processes of the higher particles, Bk ^ r Bi + s B2, where r and s are 
all those integers which satisfy 



rUk > r nil + sm2 
can be computed by Fermi golden rule 



r + s = n 



1 " 

dV = {2iTf5\P-pi p„) |T^. r 



2E J-l {2Ti)2Ei 
1=1 ^ ' 



(3.45) 



(3.46) 
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Figure 6: Neutral bound states of ip^ theory for g < 1. The lowest two lines are the 
stable particles whereas the higher lines are the resonances. 



In this formula, P denote the 2-momentum of the decay particle whereas the am- 
plitude Tfi is given by the matrix element 

Tfi = {BkiP) I . . . Bi{ps) B2{ps+i) ■ ■ ■ B2{pn) ) . (3.47) 

At the moment it is difficult to have control of this matrix element, so the best 
we can do is to estimate some universal ratios of the decay rates by phase space 
alone, assuming that the matrix elements are of the same order of magnitude for 
the various processes. For instance, taking ,^ < |, we can estimate the decay rates 

^ El + El 

E,^E, + E2 (3.48) 

^22 E^ ^ E2 + E2 

with respect to the decay rate Ff^ of the process E3 ^ Ei + Ei. We have 



^ ^ £11 ^ /m|^-4m| 



5 = £12 ^ / mj {ml - 4ml) ^ ^3 49) 

rfi V [iTT'l — (^1 — m2)'^][ml — (mi + m2)^] ' 

^ ^ ^ ma / m| - 4m^ 

rfi m^ V ml — 4m2 

Notice that at ^ = C,ci — 0.1558... the 5-particle has its mass exactly equal to 2m2. 
This makes the amplitude ratio C divergent at this point. For ^ > the decay 
process E^ ^ E2 + E2 is, instead, obviously forbidden. The plot of the quantities 
A, E and C in the range < ^ < | (for the first two) and for < ^ < ioi the 
latter, is shown in Figure d 
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Figure 7: (a): plot of decay ratios A (lower curve) and B (middle curve) as a 
function of ^ E [0, |]. (h): plot of the decay ratio C for ^ G [0,^c)- 

Since T1 = Yla t^ab is the inverse of the hfe-time of the particle Be, from 
Figure IHla one can see that the higher particle S5 is only slightly more stable than 
the lower particle B3 for C, < Cc- It is only around the critical value C,c that there is an 
enhancement of C and, correspondingly, the life-time of B^ becomes much smaller 
than the one of -B3. For C, > ^c, the ratio of the total life-time of the particles -B5 
and -B3 is given only by {A + B)^^. Plotting this quantity as a function of C,, one 
discoveres that there is a narrow window of values of given by the interval [■Ccu ■^£2]) 
with — 0.1612.., where the life-time of B^ is larger than the one of B3 (Figura 
IHlb). This counter-intuitive behavior of the life-time ratios is a simple consequence of 
the peculiar properties of the phase space in two dimensions. However, in the decay 
processes analysed in |7j, where the transition amplitudes were exactly computed, 
the tendency of the higher particles to be more stable than the lower ones was shown 
to be further enhanced by the dynamics. It would be interesting to study whether 
this is also the case of v?^ theory, using, perhaps, the numerical methods introduced 
in m. 

Let us now comment the general scenario emerging from the semiclassical anal- 
ysis. One could be obviously suspicious about the conclusion of the absence of the 
bound states for A > Ac, with the critical value Ac given in eq. ()H.42|1 : after all, 
this value is not infinitesimal and it might be, in fact, very well out of the realm of 
validity of the semiclassical approximation. This is a legitimate suspect and, in the 
absence of an exact solution which either confirms or disproves it, it is fair to say 
that the spectrum of the theory in such a strong coupling regime essentially remains 
an open question. Nevertheless, we would like to draw the attention on the following 
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Figure 8: (a): plot of the ratio of the life-time of the particles and for ^ < ^ci- 
(h): plot of the quantity [A + B)~^ in the interval C, G [0.14,0.2], with the window 
of the values of ^ where B^ is more stable than B-^ 

arguments in favour of the semiclassical analysis: 

1. Concerning Ac, in other models the critical value of the coupling is definitely 
not small either. In Sine-Gordon for instance, the critical value is given by 

= 47r. 

2. The possibility of having bound states in this theory is related to the cubic 
interaction into the potential ()3.31|) . This provides, in fact, an effective attrac- 
tive interaction between the particles which, for low values of A/m^, overcomes 
the repulsive interaction given by the v?^ term. However, the coefficient in front 
of (f^ scales as -\/A, whereas the one in front of (f^ as A. Therefore, it seems 
plausible that there should exist a sufficient large value of A/m^ where the 
repulsive force prevails on the attractive one. 

3. There is an additional little piece of information encoded in the semiclassical 
formula of the poles ()3.40|) which can be further exploited. Namely, suppose 
that the adimensional coupling constant g is very small, so that both the 
formula of the poles and the corresponding one of the masses ()3.43|) can be 
trusted. These expressions show that, by increasing g, the common tendency of 
all particles is to move forward the threshold of the two kink state and to decay 
afterward. If nothing will occur to stop this motion of the poles by increasing 
the coupling constant, the actual scenario of the theory will be then the one 
predicted by the semiclassical analysis, i.e. the existence of a critical value gc 
(surely different from the one given by the semiclassical approximation) but 
with no neutral particles beyond it. 
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4. Finally, there is an exact mapping (which will be discussed in Section be- 
tween the kinks of <y9^ and those of the Sine-Gordon model. Hence, one could 
argue that, as there is a critical value of the coupling in Sine-Gordon model in 
order to have bound states, the same should also happen for ip^ theory. 



3.1 Simple but useful observations 

In this section we will discuss some general features of the semiclassical methods 
which will be useful in the study of other models. 

We will firstly present an equivalent way to derive the Fourier transform of the 
kink solution. To simplify the notation, let's get rid of all possible constants and 
consider the Fourier transform of the derivative of the kink solution, expressed as 

/OO -1 
dxe'^^ ^ . (3.50) 
.OO cosh X 

We split the integral in two terms 

dxe"''' ^+ / dxe"'"' ^ , (3.51) 

.OO cosh X Jo cosh X 

and we use the following series expansion of the integrand, valid on the entire real 
axis (except the origin) 



cosh X 



^ = 4 y'(-l)"+ine-2'^l"l . (3.52) 

n=l 

Substituting this expression into (j3.51|) and computing each integral, we have 



n=l 



1 1 

+ 



k-2n k + 2n 



(3.53) 



Obviously it coincides with the exact result, G{k) = irk/ sinh ^k, but this derivation 
permits to easily interpret the physical origin of each pole. In fact, changing k to 
the original variable in the crossed channel, k —>■ [m — 6)/C,, we see that the poles 
which determine the bound states at the vacuum | a) are only those relative to 
the exponential behaviour of the kink solution at x ^ — oo. This is precisely the 
point where the classical kink solution takes values on the vacuum \ a). In the case 
of </9^, the kink and the antikink are the same function (up to a minus sign) and 
therefore they have the same exponential approach at x = — cx) at both vacua | ±1). 
Mathematically speaking, this is the reason for the coincidence of the bound state 
spectrum on each of them: this does not necessarily happens in other cases, as we 
will see in the next section, for instance. 
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The second comment concerns the behavior of the kink solution near the minima 
of the potentiaL In the case of <y9^, expressing the kink solution as 

m — 1 



mx) = —p= tanh 



mx 



V2 

and expanding around x = — oo, we have 



V\ ev^^ + 1 



(3.54) 



777 

ifit) = --j= [1 -2t + 2t2 -2t=^ + ■■■2(-l)"r---] , (3.55) 

where t = exp[v^x]. Hence, all the sub- leading terms are exponential factors, with 
exponents which are multiple of the first one. Is this a general feature of the kink 
solutions of any theory? The answer is positive. To prove it, consider the equation 

^ = V2t7M , (3.56) 

in the limit in which x — > — oo and (p^ifa (the same reasoning can be done, as 
well, in the other limit x +oo and ip (pf,). To simplify the expression, let's 
make the shift ?7(x) = v?(x) — (fa- Assuming that the potential U{(p) is regular near 
ifa, we can expand the right hand side of ()3.5(jj) in power of t] 

^ = aiT] + a2?7^ + a^rj^ H (3.57) 

ax 

It is now easy to see that the nature of the solution strongly depends on the presence 
or on the absence of the first term, which express the square root of the curvature of 
the potential U{ip) at ipa- In fact, if ai = a; 7^ 0, there is an exponential approach to 
the minimum, with all sub-leading terms given by multiples of the same exponential. 
To show this, let's introduce t = e^^ and express rj in power series of t, 77 (t) = 
Yl^=i f^nt"'- Substituting this expression into (j3.57|) . we have the following recursive 
equations for the coefficients Hn 

00 00 / n \ 

n=l 71=1 n=l \k=l / 



00 



E ^^k,^^k,^^k,)t'' + ■■■ (3.58) 

71=1 \ki;ki+k2+k3=n 



which iteratively permit to determine uniquely all of them. Explicitly, with the 
normalization /^i = 1, we have 

2/X2 = /i2 H /ii (3.59) 

UJ 

a2 eta 3 
3/i3 = At3 + 2 — /ii /i2 H fii 

UJ UJ 
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By summoning to well known theorems of uniqueness of the solution of the differ- 
ential equation (|3.56p . the determination of these coefficients uniquely defines the 
kink configuration near the minimum. 

When u = 0, the approach to the minimum is instead no longer exponential but 
is developped through a power-law. For instance, if the first non zero coefficient 
is ^2, by posing x = — l/(a2t), r]{t) = 7n^" and substituting into ()3.57|) . we 

have the recursive equations for the coefficients 7^ 



00 / n 



n=l n=l \k=l / 

E 7..7..7J^'^ + --- (3.60) 

n=l \ki]ki+k2+k3=n J 

which, as before, iteratively permit to fix all of them. 

The fact that the approach to the minimum of the kink solutions is always 
through multiples of the same exponential (when the curvature lo at the minimum 
is different from zero) implies that the Fourier transform of the kink solution has 
poles regularly spaced by = ^j^" ^ in the variable 6. If the first of them is within 
the physical strip, the semiclassical mass spectrum derived from the formula ()2.12|1 
near the vacuum | a ) has therefore the universal form 

m„ = 2Mab sin (n — —] . (3.61) 



As we have previously discussed, this means that, according to the value of ^a, 
we can have only the following situations at the vacuum | a): (a) no bound state 
if > 1; (b) one particle if | < .^a < 1; (c) two particles if | < < |; (d) 
■j- particles if .^a < |, although only the first two are stable, the others being 
resonances. So, semiclassically, each vacuum of the theory cannot have more than 
two stable particles above it. Viceversa, if u; = 0, there are no poles in the Fourier 
transform of the kink and therefore there are no neutral particles near the vacuum 
\ a). 

Finally, we would like to comment on the mass formula of the first neutral state, 
given by 

1 / m \ 2 



mi = 2M sin ( tt^ ) ~ m 



- 



24 \MJ 



(3.62) 



The first term coincides with the curvature of the minimum. Is there a way to un- 
derstand the presence of the subleading correction? Although it is well known (see, 
for instance, the discussion in ^H]) that this term cannot be correctly reproduced by 
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semi classical perturbation theory^, it is worth showing a simple calculation which 
indicates its relation with the dynamics of the kinks. To this aim, let's assume that 
the propagator of the kink (and the anti-kink) can be written as 

m - , (3.63) 

while for the propagator of the neutral particle we take 

^oik) - , (3.64) 

where is the curvature of the potential at the minimum. If the kink-antikink 
has Bi as bound state, there exists a non-zero 3 particle coupling and the 
possibility of a virtual process as the one shown in Figure , in which the neutral 
particle Bi splits in a couple of kink-antikink and recombines afterward. 



K 



B 



B 



K 



(-iZ) 



Figure 9: Self-energy due to the loop of the kink-antikink. 

The above diagram gives rise to the self-energy of the neutral particle. If p is its 
momentum, its explicit expression is given by 

(-.S(p2)) = (zc],^^)V(/) , (3.65) 

where ^ 

-^^P^ = j (^ q^-M^ + %e {p-q^-M^ + ie ' ^^'^^^ 

We are interested in evaluating this expression on mass-shell, i.e. at p^ = im? . The 
two denominators in J(p^) can be combined as 

1 1 /"^ 1 

dx 



q2 _M'^ + ie {p - qf - + ie Jq [g2 _ 2xq ■ p + xp^ - + ie 



l2 



^The reason is related to the intrinsic ambiguity of the normal-ordering procedure, which can 
arbitrarly alter the value of the mass scale. 
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Making the change of variable / = q — xp and introducing A(x) = —x{l—x)m'f + M 



2 



we have 



f d'H 

J{m^) = I dx 



(27r)2 (/2 - A + ze)2 

% \ % 1 Tfl 

= — / , , , = = arctan , = , (3.67) 

4:71 Jo A(x) 71 m V4M2 - V4M2 - 

where the factor i comes from the analytical continuation to the euclidean variables 
in the integral on /. In the limit .^^ ^ 0, the mass of the kink is much larger than 
the mass of the neutral particle, so that 

J(m^) ^ ^ ^ , (3.68) 



and therefore, for the value of the self-energy on mass- shell, we have 

1 

47rM2 



.S(m^)) ^.(c^_^)^-i- . (3.69) 



+ • + 



Figure 10: Propagator of the scalar particle in the ladder approximation. 

In the ladder approximation shown in Figure ^1 the propagator of the neutral 
particle becomes 

G{k) ~ Go(fc) + Go(fc)(-iS(m2))Go(fc) (3.70) 
+G'o(A;)(-zS(m?)) Go{k) (-zS(m?)) Go{k) 



^2 _ j^2 _ S(m2) 



and, correspondingly, the mass of the particle changes as 



m'^ml=m' + S(m?) = - [c'^^^f^l^ . (3.71) 

Notice that this correction is always negative, i.e. the presence of the kink tends 
to decrease the value of the mass of the neutral particle, initially expressed by the 
curvature of the potential. If it happens that, by varying the coupling constant, 
the second term exceeds the first, there is an imaginary value in the mass mi. 
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This implies that the particle disappears from the stable part of the spectrum. 
Considering the coefficient c]^^ as fixed, this occurs for sufficiently large values of 
the coupling constant. Viceversa, decreasing the coupling, the mass of the kink 
becomes larger and the correction gets consequently smaller, making the mass of 
the neutral particle closer to its perturbative value. 

3.2 Watson's equation and the ^-matrix 

As a final topic of this section, we now consider the issue of the Watson's equation 
satisfied by the Form Factor. This concerns, in particular, the interesting possibility 
of extracting the S'-matrix of the kinks, at least in certain regimes of rapidity and 
coupling constant. The basic idea behind the Watson's equation is simply the com- 
pleteness of the asymptotic states: given a matrix element of a local operator Q on 
a given in state | i.e. F£ = (0 | ^(0) \n )in, one can employ the completeness 
relation of the out states to get the following relation^ 

Fin^ = (0 I ^ I n)in = ^(0 I g I m)outout{m \ n)in 

m 

m m 

In this equation Sn^m is the S'-matrix amplitude relative to the scattering of the n 
initial particles into m. Let's suppose now that the momenta of the incoming par- 
ticles are so small that it is impossible to open higher inelastic channels. Moreover, 
let's assume that are also absent "decay" processes in lower mass particle states. If 
this kinematical regime exists, then the remaining non-zero scattering amplitudes 
in fl3.72|) are nothing else but elastic, so that the Watson's equations become similar 
to those employed in the integrable models fH], i-e. 

(0 I G{Q) I n),„ ~ Sn^n (0 I G{Q) I n)out . (3.73) 

In this case, from the ratio of the F^^V-^i"]' could get the elastic part of the 

S'-matrix S'„_»„, an expression obviously valid only in the kinematical region below 

the lowest threshold of the ra-particle channel. Concerning the two-body S'-matrix, 

it is important to notice that its elastic part is a pure phase only for real values of 

9 below threshold. In the following, however, we enforce it to be a phase also for 

complex values of the rapidity since, anyhow, this is the best we can do to obtain 

an estimate of this quantity. 

^This is a scheleton form of the Watson's equation: the matrix elements depend on the momenta 
of the particles and therefore the sum of the intermediate states stays also for a multiple integral 
on these variables. 
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Let's follow this suggestion to see whether it would be possible to determine the 
elastic part of the S matrix of the two kink states in the (f'^ theory. First of all, we 
have to establish that there are no neutral particles i?„ with mass m„ < M otherwise, 
for the non integrability of the theory, the "decay" channel | KK) ^\ BnBn) will 
always be open, even if the kinks are at rest (here K denotes the anti-kink). Since 
rrin > TTLi, it is sufficient to impose mi > M in order to prevent such decays. This 
gives rise to the following condition on the coupling constant 

sin^>i z.e. e>^- (3.74) 

Once we are in this range of the coupling constant, the absence of the higher mass 
thresholds is ensured by taking sufficiently small values of the rapidity difference of 
the two kinks. 

Notice that the Watson's equations are valid irrespectively of the operator Q, 
an important point on which we shall come back later. Taking for granted this 
insensitivity to the operator Q, then we can take the Form Factor of the field (p{x) 
computed in (j3.39|) . The Form Factor of the out state is simply obtained by substi- 
tuting in ()3.39j) 9 — > —9. By the ratio of these quantities, we arrive to the putative 
expression 

sinh ^ 



Slt-aiO) = S-r''{9) ^ S{9) . (3.75) 

sinh i^^^j 

If correct, this formula should describe the elastic scattering of the kinks 

below their inelastic threshold. 

The two amplitudes of the kink scattering can be represented by the diagrams 
of Figure ^2 where the indices on the left and on the right are those relative to 
the initial and final vacua respectively, whereas the indices on the top and on the 
bottom are the other vacua "visited" during the scattering^. 

The amplitudes obtained in eq. (|3.75p satisfy (by construction) the unitarity 
equation 

S{9)S{-9) = 1 . (3.76) 

However, they fail to satisfy the crossing relation which is expected for the correct 
S'-matrix 

SZ-a i'^ - ^) = S~:'-^{9) . (3.77) 

®In this particular theory, it is impossible to change the intermediate vacua in the scattering 
process: this means that the kinks of (^^ essentially behave as ordinary particles. 
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Amplitude 



s-channel 



t-channel 



Figure 11: Elastic scattering amplitudes of the kinks and their poles in the s-channel 
and in the t-channel. The dots at the vertices are the on-shell 3-particle couplings. 

The validity of this relation simply follows by turning of 90*^ one of the amplitudes 
of Figure ^2 and comparing with the other. In addition to this problem, the above 
expression for S has another drawback, as it becomes evident by a closer look at 
its analytic structure: even though it has all the poles of the bound states in the 
s-channel correctly localised at 6* = Z7r(l — n^), their residue 



is not always (imaginary) positive as, instead, it should be. In the vicinity of these 
poles the correct S-matrix should indeed reduce to 



and, for unitary theories, the 3-particle coupling (jf"^ ^ is real. 

To make a further progress toward the correct identification of the elastic S ma- 
trix of the kinks, let's now explore the arbitrariness of the operator Q{x) entering the 
Watson's equation. Semiclassically, for any operator which is a reasonable function 
G{(f) of the field (f{x), its Form Factor is given by 



fab = {KaM I G[ip{0)] I K^,ie2) = / dxe^^-^^Gb„,(x)] . (3.80) 




(3.78) 



S{9) = i 



- iu\a 



(3.79) 




oo 
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To be defined, let's consider a class of operators expressed by power series in 



Q = ai^p + a2^ H ctn V^" H (3.81) 

The regular part of the above Form Factor is obtained in terms of the derivative of 
the function inside the integral which, in a simplified notation, is given by 



oo 



ikx dG fd^\ 



oo 



H(k)= I d.e'"- ^-^j , (3.82) 



where ip denotes here the kink solution. It is easy to see that H{k) has always the 
same poles of the Fourier transform of (^) (just using the previous argument on 
the asympotic behavior of the integrand). However, its residues can be arbitrarirly 
varied by changing the coefficients a„ of the expansion ()3.81|) . Said in another way, 
at the semiclassical level, the most general expression of the Form Factors for the 
(yj^ theory is given by 

oo 
n=— oo 

where the coefficients s„ are arbitrary numbers. By varying them, the only effect 
is to change the position of the zeros of the function H{k), obviously leaving the 
position of its poles untouched. 

Reestablishing now the original variable k {in — 9)/S, and taking the ratio of 
the Form Factors of an operator of the above class 

we see that the resulting function can be any function which fulfills the unitarity 
equation 

R{e)R{-e) = 1 , (3.85) 

and which has the correct poles in the s-channel. Arbitrary other poles of R{0) 
may come from the zeros of if (ztt + 6). In conclusion, the farthest we can go in 
the use of the semiclassical form factors, is to fix the poles of the S-matrix in the 
s-channel alone. To find its actual expression, one necessarily needs to integrate 
this information with the others coming from the physical nature of the problem at 
hand. 

For the ip'^ theory this is easy. In fact, an S'-matrix which must simultaneously 
satisfy the two equations 

*(«^(-« = ' ■ (3,86) 

S{P) = S{i7r -9) , 
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can only be expressed as a product of the elementary functions 

tanh h(6 + mri) 
[rj] = [l-rj] = -4)- — ^ . 3.87 

These functions have two poles, with residues of opposite sign: one at ^ = ivr?], the 
other at its crossing symmetric position 9 = m{l — rj). Hence, the natural proposal 
for the S'-matrix of the kinks below their inelastic threshold is 

n 

SZM = (-l)"+^n[^^] ' (3.88) 

k=l 

an expression which holds for the following values of the coupling constant 

^ <i<- ■ (3.89) 
n + 1 n 

By construction, it satisfies both the unitarity and the crossing symmetry equations. 
It has the right poles in correspondence with those expected from Figure El It is 
also easy to check that it has the correct positive residue at all poles in the s-channel 
(and a negative one at all the t-channel poles) . 

In view of the condition ()3.74|) . the only physical values that n can take in the 
above expression are n = 0,1,2. Somehow surprisingly, it keeps its validity (for 
instance, the positivity of the residues at all poles) even for other values of ^, where 
there are decay processes into lighter breathers. Notice that, when there are no 
breathers in the theory (n = 0), the S'-matrix of the kinks simply coincides with the 
one of the Ising model in its low temperature phase, S = —1 



3.3 Finite volume 

A way to investigate the spectrum of a quantum field theory is by studying its 
euclidean version on a finite volume, say on an cylinder of width R along the space 
direction and infinitely long in the euclidean time direction r. Assuming periodic 
boundary conditions ip{0,T) = (p{R, r), how the spectrum Ei{R) of the finite- volume 
Hamiltonian would look like? 

For theory, the answer is as follows. First of all, for dimensional reasons, the 
energy levels can be cast in the form 

271 

Ei{R) = —e,{MR) , (3.90) 
rC 

where ei{MR) are the scaling functions'^ . They depends on the adimensional pa- 
rameter MR, where M(A) is the mass of the kink, which we assume to be always 

^The semiclassical expression of the scaling functions of ip^ theory with anti-periodic boundary 
conditions has been studied in |14j . An issue of a particular interest is their crossover from the 
conformal regime to the massive behavior. 
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finite. Secondly, it is easy to foresee tlieir beliavior in two limits, i? — and R — >■ oo. 

For i? — ^ 0, the theory presents a conformal invariance. The scaling functions 
becomes then ej(0) = (2Ai — c/12), where c = 1 is the central charge of the bosonic 
theory, whereas Aj are the conformal dimensions of the various conformal fields 
present in this limit. In the other case, R oo, the theory displays instead its 
massive behavior. Therefore all levels go to a multi-particle state, with a mass gap 
A4 given by the sum of the masses of the excitations entering this state. They can be 
either breathers or kinks (for periodic boundary conditions, there are only couples 
of kink-antikink states). Taking into account a possibile bulk energy term eo(A), the 
energy levels are then expected to go asymptotically as 

Ei{R) eo R + Mi , R-^oo (3.91) 




0.5 1 1.5 2 2.5 3 3.5 



Figure 12: Energy levels of the bound states, as functions of MR, for the symmetric 
well potential in a finite volume with periodic boundary conditions. 

At a finite volume, however, there is a finite energy barrier (order R) between the 
two vacua and they will be in contact each other through the tunneling of the kink 
states. Correspondingly, the energies of the vacua, together with all the energies of 
the excitations above them, have an asymptotical exponential splitting: they come 
in pairs and become doubly degenerate only in the infinite volume limit. A typical 
outcoming of this circumstance is shown in FiguralT^ 



27 



To check whether the the semiclassical prediction is correct or not, it would 
be sufficient to study the movement of the energy levels by increasing A. If the 
prediction of a critical value beyond which there are non longer bound states is 
correct, one should observe a progressive approach of all their energy lines toward 
the two-kink threshold 2M and their disappearance into the continuum once A > Ac- 



4 Asymmetric wells 



In order to have a polynomial potential with two asymmetric wells, one must nec- 
essarily employ higher powers than ip^. The simplest example of such a potential is 
obtained with a polynomial of maximum power (p^, and this is the example discussed 
here. Apart from its simplicity, the (p^ theory is relevant for the class of universality 
of the Tricritical Ising Model. As we can see, the information available on this model 
will turn out to be a nice conffimation of the semiclassical scenario. . 

A class of potentials which may present two asymmetric wells is given by 

U{v) = -[^ + a^\ U - 6^ [v^ + c—] , (4.92) 



2 V ' ^VxJ V WxJ V A 

with a, b, c all positive numbers. To simplify the notation, it is convenient to use the 
dimensionless quantities obtained by rescaling the coordinate as — > mx^ and the 
field as (p{x) —>■ VX/mip{x). In this way the lagrangian of the model becomes 



r = — 

A2 



(4.93) 



The minima of this potential are localised at ip^Q^ = —a and ipf*^ = b and the 
corresponding ground states will be denoted by | ) and | 1 ) . The curvature of the 
potential at these points is given by 

U"{-a) ^ uj', = {a + bna' + c) ; 

U"{b) = Lul = (a + 6)2(62 + c) . ^ ' ' 

For a ^ b, we have two asymmetric wells, as shown in Figure El To be definite, 
let's assume that the curvature at the vacuum | ) is higher than the one at the 
vacuum | 1 ), i.e. a > b. 

The problem we would like to examine is whether the spectrum of the neutral 
particles \ B)^ {s = 0, 1) may be different at the two vacua, in particular, whether it 
would be possible that one of them (say | 0)) has no neutral excitations, whereas the 
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Figure 13: Example of Lp^ potential with two asymmetric wells and a bound state 
only on one of them. 

other has just one neutral particle. The ordinary perturbation theory shows that 
both vacua has neutral excitations, although with different value of their mass: 



= (n + b)^/2{^i^Tc) , m(^) = {a + b)^/2 (&2 + c) . (4.95) 



Let's see, instead, what is the semiclassical scenario. The kink equation is given 
in this case by 

^ = ±{^ + a){ip - b) ^^^ + c . (4.96) 
ax 

We will not attempt to solve exactly this equation but we can present nevertheless its 
main features. The kink solution interpolates between the values —a (at x = — oo) 
and 6 (at X = +oo). The anti-kink solution does viceversa, but with an important 
difference: its behaviour at x = — oo is different from the one of the kink. As a matter 
of fact, the behaviour at x = — oo of the kink is always equal to the behaviour at 
X = +00 of the anti-kink (and viceversa), but the two vacua are approached, in this 
theory, differently. This is explicitly shown in Figure IT^ and proved in the following. 

Let us consider the limit x — — oo of the kink solution. For these large values 
of X, we can approximate eq. ()4.9(i|) by substituting, in the second and in the third 
term of the right-hand side, ip ~ —a, with the result 

(^j ^ (v9 + a){a + b)yja^ + c , x ^ -oo (4.97) 
\ax J 

This gives rise to the following exponential approach to the vacuum | 0) 

(pQ^i[x) ~ —a + y4exp(co'oa;) 7 x — > —00 (4.98) 

where A > is a arbitrary costant (its actual value can be fixed by properly solving 
the non- linear differential equation). To extract the behavior at x ^ —00 of the 
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Figure 14: Typical shape of (^)q^, obtained by a numerical solution of eq. \4.9(^ . 

anti-kink, we substitute this time cp c:^ b into the first and third term of the right 
hand side of ()4.96p . so that 

dip 



~ {ip-b){a + b)Vb'^ + c , x->-oo (4.99) 

/ 1 

This ends up in the following exponential approach to the vacuum | 1) 

ipifi{x) b — B exp{uJix) , X ^ —oo (4.100) 

where S > is another constant. Since c^o 7^ u^i, the asymptotic behaviour of the 
two solutions gives rise to the following poles in their Fourier transform 



ujq + ik 

(4.101) 

^{^1,0) ^ 

In order to locate the pole in 6, we shall reintroduce the correct units. Assuming to 
have solved the differential equation ()4.96p . the integral of its energy density gives 
the common mass of the kink and the anti-kink. In terms of the constants in front 
of the Lagrangian ()4.9H|1 . its value is given by 

171 

M = —a , (4.102) 

where a is a number (typically of order 1), coming from the integral of the adimen- 
sional energy density ()1.5|1 . Hence, the first pole^° of the Fourier transform of the 

^°In order to determine the others, one should look for the subleading exponential terms of the 
solutions. 
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kink and the antikink solution are localised at 



m \ , f A 



iTi[l-u}o — - ] = m [l-ujQ 
V vrM/ V am 



,4 



(4.103) 



If we now choose the coupling constant in the range 

1 1 

- < ^ < - , (4.104) 

the first pole will be out of the physical sheet whereas the second will still remain 
inside it! Hence, the theory will have only one neutral bound state, localised at 
the vacuum | 1 ) . This result may be expressed by saying that the appearance of a 
bound state depends on the order in which the topological excitations are arranged: 
an antikink-kink configuration gives rise to a bound state whereas a kink-antikink 
does not. 

Finally, notice that the value of the adimensional coupling constant can be chosen 
so that the mass of the bound state around the vacuum | 1 ) becomes equal to mass 
of the kink. This happens when 

4 = ^ ■ (4.105) 

Strange as it may appear, the semiclassical scenario is well confirmed by an 
explicit example. This is provided by the exact scattering theory of the Tricritical 
Ising Model perturbed by its sub-leading magnetization. Firstly discovered through 
a numerical analysis of the spectrum of this model [25, its scattering theory has 
been discussed later in |22]- It involves several amplitudes but, for our purposes, 
it is enough to focus only on those given below. With the same meaning of the 
diagrams as in the previous section, their exact expression is given by 
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The function So{0) ensures the unitarity condition of the whole set of amphtudes 
and it is given by 

-I) ^ To) ^ To) t i) t -1) t I) t {0, -I) 



w 



sinh^( 



sinh|^(^^-^) 



(4.106) 



where 



w{9, x) 



sinh {^0 + inx) 
sinh (^6^ — zTTx) 



sinh i 

t(e,x)- 2 



sinh i( 



The structure of poles and zeros of the S-matrix of this problem is quite rich but, 
on the physical sheet, <lm 9 < in, the only poles of the S'-matrix are located at 
^ = ^ and ^ = f . 





00 
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s-channel 
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Figure 15: Elastic scattering amplitudes of the kinks in an asymmetric wells potential 
and their intermediate states in the s-channel and in the t-channel. 

The first pole corresponds to a bound state in the s-channel whereas the second one 
is the singularity due to the particle exchanged in the crossed t-channel. However, 
the residues of the two amplitudes at 6* = ^ are quite different! In fact, for the 
first amplitude we have 



Resg^2j^ 'S'oo(^) = ; 



while for the second 



ReSfl 



2iTi 
■ 3 



ill 



(4.107) 



(4.108) 
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where 



UJ 



5 ^ (I) ^ (l^) ^ (I) ^ (I) (^) 



(4.109) 



= sin(7rx)). Hence, in the s-channel of the amphtude 6*00, there is no bound 
state related to the vacuum | ) : its only singularity comes from the bound state 
on the vacuum | 1 ), exchanged in the t-channel. In the amplitude the situation 
is reverted (the two amplitudes are related by crossing): there is the s-channel 
singularity due to the bound state present on the vacuum | 1 ) while the one of the 
t-channel is absent. This is easily seen in Figure ITH| where the original amplitudes 
are streched along the vertical direction (s-channel) and along the horizontal one 
(t-channel). 

A simple way to check the above scenario would be to study the finite volume 
energy spectrum with periodic boundary conditions for the field (y9(x, r) at the edge 
of the cylinder of width i?, as was done, in fact, for the Tricritical Ising Model in 
PT] . At a finite R the energies of the two vacua are exponentially splitted through 
the tunnelling process of the kinks. However, this does not occur for their excitation. 
If we consider, for simplicity, the case of only one bound state, its energy level is 
a single, isolated curve placed between the vacua energy and the threshold lines 
(FiguralT^. This situation has to be contrasted with the one shown in Figura IT^ 
related to a potential with two symmetric wells. 




Figure 16: Energy levels of the asymmetric well potential in a finite volume with 
periodic boundary conditions. 
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5 Sine-Gordon model in the semiclassical limit 



In this section we will use the exact solution of the Sine-Gordon model in order to 
check the semiclassical approximation and to learn some important lessons from this 
comparison. The potential is given in this case by 

2 

171 

U{v) = _(l-cos/5^) . (5.110) 

There is an infinite number of minima, ipn'^ = 27m/, which correspond to the quan- 
tum vacua \ n). For their equivalence, one can choose to study the excitations on 
one of them, say the vacuum | ). 



5.1 Exact scattering theory 

The exact scattering theory of this model has been discussed in P and it will be 
briefly summarised below. The kink of the Sine-Gordon model interpolate between 
the vacuum | ) and its neighbouring ones | ±1 ), and their scattering processes are 
described by the amplitudes 

I Ka,o{ei)Ko,-ai02)) = | i^a,o(^2)i^O,-a(^l)) (5.111) 

with a = ±1. In the neutral kink-antikink channel, the amplitude Sqq{6) = Sr{6) 
describes their reflection process while >S'qq~''(^) = St{0) describes their transmission. 
In the kink-kink scattering there is only the transmission amplitude >S'^°_^(^) = S{6). 
The reason of this terminology stays in the identification of the states, due to the 
equivalence of the various vacua: for instance, the kink | K_i q) must be identified 
with I Kq^i) and similar identification can be also established for the others. The 
above amplitudes, represented as in Figure [13 satisfy the unitarity equations 

Sr{9) Sr{-9) + St{9) St{-9) = 1 ; 

Sr{9) STi-9) + STi9) SRi-9) = ; (5.112) 
S{9)S{-9) = 1 , 

and the crossing symmetry relations 

Sr{z7i-9) = Sr{9) ■ ST{^n-9) = S{9) . (5.113) 
Their closed solution is given by P 

sinh ^9 

St{9) = . , S{9) ■ 

smh j[in — 9) 
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TV 

^ gl]2 - — 



where ^ is the so-called renormalised coupling constant 

1 



while the amplitude S{9) is given by 



57r 1 _ 



(5.115) 



Stt 



S{9) 



exp 



—I 



dt sinh |(7r - 
t sinhf cosh^ 



sin 6t 



(5.116) 



The pole of S{6) are all outside the physical sheet, as it can be read from its equiv- 
alent infinite-product representation 



3(6) 



n 



r((- + i)i + ^l^)r(i + (n + i)£-zA) 



X 



!)r((n + i)4-.|,)r(i + (n + i)^ + ^A) 



(5.117) 



Hence, the bound states can be obtained from the poles inside the physical strip of 
the amplitudes Sr and St- Since they are located at 



k = l,2,... 



(5.118) 



^ must be less than 1 in order to have a bound state. This leads to the critical value 
of the coupling constant given by 



(5.119) 



beyond which, there are no bound states. For (3 < j3c, their number is given by 



N 



. Calling M the mass of the kink, their mass is expressed as 



rrir. 



2M sin n - 
2 



n = 1,2,...A^ 



(5.120) 



Obviously, the stability of the particles with mass m„ > 2mi is ensured in this case 
by the integrability of the model. Concerning the residues of the amplitude Sr and 
St in the s-channel, they are alternatevely equal or opposite, depending on n 



Res 5/?[z7r(l-0] = (-1)" Res 5T[«vr(l - <)] = (-1)"+^- sin ^ 5[z7r(l - <)] • 

(5.121) 

To appreciate the meaning of this result, let's introduce the combinations which 
diagonalise the S'-matrix and which display, in this case, also the charge conjugation 



TT 
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Figure 17: Elastic scattering amplitudes of the kinks and their poles in the s-channel 
and in the t-channel. 

symmetry of the kinks 



(5.122) 



, , cosh ^(ZTT + 9) 

5'->W= MO) = 3(0) 



In view of (j5.12ip . the bound states with n odd appear as poles only in S^~\9) 
whereas those with n even appear as poles only in S^'^\9). The two sets couple, 
respectively, to the combination of the kink states 

I /io,iA'o,i) ± I A'o,-iA'_i,o) . 
The above result permits to check explicitly the non-degeneracy of each mass level 
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(j5.120p . As a matter of fact, the equal or opposite value of the residues of Sr 
and St is a general feature, which is valid each time that a vacuum state is in 
communication with two symmetric neighbouring vacua. It is worth spending few 
words on its derivation and on its possible generalization. 




Figure 18: Portion of a potential with three neighbouring vacua, two of them sym- 
metrically placed with respect to the central one, which supports a certain number of 
bound states. 

5.2 Residue relations 

Consider a potential which has, locally, a situation like the one shown in Figure ITHl 
The vacuum in the middle, which we denote by | ) , exchanges kinks with the two 
neighbouring ones, here labelled by | ±1 ). The latter vacua are equal between each 
other, but not necessarily equal to the central one. This means that the kinks | Kqi) 
and I -ft'c-i) (as well as their anti-kinks) have the same masses, and that | KqiKqi) 
has the same properties of | ii'o,-i-f^-i,o)- In this situation, even if the theory may 
be not-integrable, one can define an (elastic) reflection and transmission amplitudes 
for the scattering of the above states, just as in the case of Sine-Gordon. Obviously, 
in the non-integrable case, this picture will only be valid in the kinematical region 
below their inelastic threshold (and also with the lowest mass of the bound states 
higher than the mass of the kink themselves). 

Suppose now that the vacuum | 0) has certain bound states | i?„ )o- Let's 
call (7^^ and g^^^ the coupling to the n-th bound state to the kinks | KqiKiq) and 
I KQ _iK_ifi) , respectively. The bound states show up in the poles of both Sr and 
St- Streching the corresponding amplitudes along the s-channel, the residue in St is 
proportional to {g^^^g^+^), whereas the residue in Sr is proportional either to {g+^Y 
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Figure 19: Reflection and transmission amplitudes for the kinks \ KoaKao) (a = ±1 ). 
The black dot represents g^]f^ while the white one g^^'^ . 

or to ((^i"^)^, depending on which amphtude one is looking at (Figure E)). However, 
for the equivalence of the two vacua | ±1 ), the last two quantities must be equal. 
Hence 



i.e. the residues of the and St are either equal or opposite. Correspondingly, 
at the pole, the even, or the odd combination of the amplitudes, becomes a one- 
dimensional projector: the neutral particle state is not degenerate, as we already 
know. 

The reasoning can be further generalised in the case of three neighbouring vacua, 
one different from the other, but with kinks of the same masses. This is not im- 
possible: suppose, in fact, that we deform the potential of Figure IT^ but always 
keeping equal the integrals ()1.5p of the right and left kink. Notice that, by using the 
equation ()1.3|1 . eab{x) can be equivalently expressed as 




±9- 



(n) 



(5.123) 
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So, we are simply looking for a potential which supports two kink configurations 
(po^±i{x) that, although different, satisfy however 




(5.124) 



There is, of course, no mathematical obstacles in fulfilling this request^^. Concerning 
the shape of the potential, it may look like the one shown in FiguraEUl 




Figure 20: Portion of a potential with three neighbouring vacua of different shape 
but with kinks of equal mass interpolating between them. 

Under this deformation, the masses of the bound states may change but their 
values can never cross, for the non-degeneracy of these particle. The kinks, on the 
contrary, are still degenerate and therefore they can mix each other under the scat- 
tering processes. Due to their asymmetric shape, the elastic part of their scattering 
is described, in this case, by 3 amplitudes. Using the short notation K± and K± to 
denote the left and right kink/ant ikinks, we have 

\ K+{9,)K+{92)) = S2{9) \K.{92)K.{9^)) + \ K+{92)K+{9^)) . 

(5.125) 

Calling as before (yfi"^ and g^]f^ the on-shell coupling of the left and right kinks to the 
bound states, the residue r2 of 5*2 will be proportional to (g^^^g^^^), whereas those 
of 5*1 and S'3 (ri and r^) will be proportional to (fi'i"'')^ and ((7+"'')^, respectively. 
However, the residues ri and are not equal in this case, although they satisfy the 
relation 

rirs = (ra)' . (5.126) 

"'^"'^Such a potential can be found, for instance, by empfoying the deformation procedure described 
in EH]. 
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Diagonalising the above S'-matrix, we have two different phase shifts 

^^5^ie) ^^(^g^^g^^ ^ ^s, - S,)^ + AiS^)'^ . (5.127) 

Accordingly to the relative sign of g^^^ and g^_!^\ one of the two will have a vanishing 
residue of the pole of the bound state, as it is easily checked by using eq. ()5.126j) . 
The other amplitude will act, then, as a one-dimensional projector. 

An explicit example of an asymmetric 3-vacua kink scattering is provided, for 
instance, by the following set of functions 

where 5*0 (6*) is the same of eq. ()4.106p and ensures the unitarity condition of these 
amplitudes. It is easy to see that, at ^ = 27rz/3, the residues of Si and 5*2 are 
different, although they satisfy the relation (j5.126|) . 



5.3 Semiclassical analysis 



Before proceeding with the analysis of the classical kink configurations of the Sine- 
Gordon theory, let's underline one fact that will be important later. Namely, if we 
take the limit ,^ — > in the exact scattering amplitudes, the S-matrix develops an 
essential singularity in C,- Moreover, it is no longer a meromorphic function of 6. 
The term responsable for this singular behavior is So{6), with the breaking of the 
analiticity of this function due to the accumulation of its infinite number of poles 
and zeros when ^ — 0. In this limit, its explicit expression can be obtained by using 
the infinite product representation ()5.117p . It reads as follows [H ITU I 



S{9) 



exp 



log 



-ITI 



di] 



(5.129) 



Turning now the attention to the semi-classical data of the kinks, the solutions 
of ()1.3j) are given by 

cj)a,a±i{x) = ±1 arctan(e"^^) , (5.130) 

with the usual identification of the vacua, modulo 2tt/(3. Since the kink and anti- 
kink solutions are equal (up to a sign), the semiclassical spectrum on | ) and on 
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its neighrouring vacua are the same, as it is expected from their equivalence. The 
classical mass of the soliton is given by Including its first quantum corrections, 
one has 

8m m m , , 

M = — = — . (5.131) 



m 



Let's now compute the semiclassical Form Factor of an even and an odd operator. 
For the even one, we can choose, for instance, G+ = A+ cos[Pip{x)] 



dx e 



:Mex 



TT 



A, 



M 



m J sinh ^ 



cos[/5<y5oi(a;)] 



(5.132) 



while for the odd one, G_ = A- sm[Pip{x)] 



A^ {K, 



sin[/3(^(0)] I iToi) = X 



dx e 



iMBx 



sm[p(foi{x)] 



-ITT A- 



m J cosh 



(5.133) 



2? 



In the final expressions we have isolated the term 6*72. The reason comes from a close 
comparison with the exact expression of the Form Factors of the above operators, 
which is available in the literature fHl- By choosing a proper normalization A± of 
the operators, they can be written as 



/W(^) 



9 



2 sinh ^ 



2 cosh 2| 



' ^ fcxlct (^) 



sinh 



9 



2 sinh ^ 



fi^lM = sinh 



2 cosh 2| 



(5.134) 



where 



Fmin{0) = exp 



dt sinh 1(1 -0 sin^i: 



(5.135) 



t sinh ^ cosh | sinh t 

We see that, a part of the function F^ij^{9), on which we are going to comment soon, 
the remaining expressions can be made equal by substituting 9/2 sinh6'/2, an 
operation which is definetly permitted within the semiclassical approximation. 

From a numerical point of view, as far as ^ < 1, the above substitution is really 
harmless. In fact, at fixed ^, the function Fmin(^) asympotically goes as 



-^mm(^) ^ exp 



^ (1-0 

4 e 



^ ^ oo 



(5.136) 
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and, therefore, for ^ < 1, the ratio F-^auiiG)/ sinh ^ becomes neghgible before one has 
the possibihty to appreciate the difference between the term 6/2 and sinh(^/2): plot- 
ting together the semiclassical and the exact expressions, one can hardly distinguish 
them on the entire infinite range of 9 fFigurel^Tj). 




Figure 21: Comparison at ^ = 0.3, between the exact form factor of cos j3ip{x) 
(continous line) and its semiclassical expression (dotted line). 

One may have noticed that the form factors f^^\0) of the operators cos (3{^p) 
and cos I3{ip) has a finite value at 6* = 0. The reason is that they are local fields 
with respect to the kinks. The same is true for all integer harmonics cos{n(3ip{x)) 
and sin(n/9<y9(x)). In fact, the semilocal index 7^ of the exponential operator e*"'''*^^^ 
with respect to the kink is [H] 

7„ = ^ . (5.137) 

For the form factor g°'{6) = {K{6i) \ e'^°"'^^^'^ \ K(62) of these operators, 7 is the 
quantity that rules their residue at 6* = 

Rese=o^"(^) = (1 - e^"^^") (0 | e^°^(°) | 0) . (5.138) 

However, taking the operators cos f v^(a;) and sin it is easy to see that the form 
factor of the first operator should have a pole at 6* = 0, whereas the other should not. 
This is, indeed, well confirmed by the semiclassical expression of these quantities. 
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For the first, we have in fact 



2 



(J^Ol(^^l) I COS 

M 1 



Km) 



dxe 



iM e X 



COS 



2'/'oi(a;) 



—Z TT 



m sinh^:^ 

2m 



^ sinh2| 



(5.139) 



while for the second 



M 



sm 



1^(0)' 



Koi) 
1 



dxe 



iMBx 



sm 



TT 



m cosh 



-IT Me 

2m 



^ cosh^i 



2<^oi(a;) 



(5.140) 

m — 6*, it is 
with 



Concerning the bound states, making the analytic continuation 9 - 
easy to see that one obtains the exact spectrum of the bound states 
n even, from the poles of the form factors of the even operators, as for instance 
F^^\6) = f^^^iir — 9),. Viceversa, from the poles of the odd operator form factors, 
like F^^\6) = /*■"'' (ivr — 9), one obtains the exact mass of the particles with n 
odd. This happens because the neutral bound states are eigenvectors of the charge 
conjugation of the model, with eigenvalues (—1)" and, therefore, the two sets couple 
only to operators with the same parity. 

Let's now discuss the issue of the 5- matrix. Taking, for instance, the ratios of 
the semiclassical form factors of the operators cos ^(f and sin |<yC computed iir =F ^5 
one obtains 



Si-) 



(+), 



ITT — 



2 

2 

F[~\i7r + 



sinh ^(«7r + 
sinh ^ {in — 



(5.141) 



cosh ^{in 



0) 



cosh ^ [in 



0) 



Comparing these results with eq. ()5.122|) . we see that they are remarkably close to 
the exact expressions but, nevertheless, they miss the important term S{(3). This 
is related to the problem, previously discussed, of the impossibility of recovering 
the correct crossing symmetry of the 5* by using the semiclassical Form Factors of 
a generic operator. In this case, however, this problem seems to get even worse for 
the mathematical nature of the function which is missing, eq. ()5.129p . In fact, this is 
a function with a branch cut in 9, that can only be obtained by an accumulation of 
infinite number of poles and zeros. Together with its essential singularity at ,^ = 0, 
its non-analytic behavior makes S{9) completly invisible to the semiclassical formula 
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(|2.12|) we are using. Saying differently, the term S{6) could only be recovered by 
the ratio of the functions F^i^{i7r±6). By the exact solution of theory we have 
indeed 



(in + 



S{e)F^S'^-e) , F^UO) = i^mi„(-^) . (5.142) 

However, in order to obtain the expression ()5.129|) . such function should behave, for 
C ^ 0, as 

2- 



Frain{0) = exp 



k=l 



1 



7r(2fc+2-r?) 



(5.143) 



dr] log 

1 ^ / 

^7r(2fe+2+r/) 

It is easy to convince oneself that the semiclassical formula ()2.12|) can never produce 
such an expression for the form factor of any operator G{ip{x)) which is an analytic 
function of the field (f{x). 

In the end, we would like to close this section with a remark that may lighten 
the bound states of (f^ theory, analysed in Section IHl To do so, it is convenient to 
rescale the coordinates x'^ — >■ mx'^ and the field (f — f3(f of the Sine-Gordon theory, 
so that its Lagrangian becomes 



r - — 

/?2 



-{d<^f + (cos, 



Under the substitution 

(p±{x,t) = TT ± 4 arctan $(a;, t) 
the Lagrangian of the Sine-Gordon model becomes 



2\2 



(5.144) 



(5.145) 



(5.146) 



The static solutions of this Lagrangian coincide with those of ^9^, once the coupling 
constants of the two theories are related as 

A P 



(5.147) 



By inserting in this relation the critical value of the coupling of the Sine-Gordon, 
/3c = v^47r, we may get a different estimate of the critical value of ^p^ theory 

Ac 



^ = Vvr = 1.77245. 



(5.148) 



This value is of the same order of magnitude of ()3.42|) . previously obtained. It would 
have been, probably, too ambitious to search a better agreement between the two 
values. After all, even though the static solutions of the two models are similar, 
their time dependent solutions are instead rather different (the Sine-Gordon has an 
integrable dynamics, whereas ip^ does not). 
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6 Double Sine-Gordon wells done doubly well 



In this section we consider the situation of a vacuum state in communication with 
two neighboring ones through kinks of different masses. A prototype of this situation 
is given by the Double Sine-Gordon model, with potential given by 

V{^) = -Acos/3(^-Acos(^|^ + 5^+C. (6.149) 

By choosing 

and by varying A, the shape of the potential changes as shown in Figure 





Figure 22: Shapes of the Double Sine Gordon potential by varying A. 

There are two regions, qualitatively different, in the space of parameters, the 
first given by < A < 4// and the second given by A > 4/i. They are separated 
by the value A = 4/i (where the curvature at the minima vanishes) which can be 
identified phase transition point [S]. 

Let's focus our analysis in the coupling constant region where A < 4fi. Switching 
on A, there are several effects on the potential: the original minima of the Sine- 
Gordon, located at ipmm = 0; ^ (i^od^), remain degenerate but they move to 
^min = — </2o! ^ + V^o (ii^od^), with LfQ = | arcsiu The different minima has, 
however, the same curvature, given by 

1 A^ 

= /i . (6.151) 

16 fi 

In correspondance with the above shifts, there are two different types of kinks, one 
called "large kink" and interpolating through the higher barrier between —00 and 
^ + 00, the other called "small kink" and interpolating through the lower barrier 
between —^ + (f>o and —(po- Their classical expressions were explicitly given in the 
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reference 



= ^ + ^ arctan 



, , 71 A 

^s{x) = + ^arctan 



'4/i + A , fm 

tann — x 

4/i - A V 2 



47r 

(mod — ) , (6.152) 
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'4/i — A /m 

tann — x 

4/i + A V 2 . 



47r 

(mod — ) 



(6.153) 



For the following considerations, we can neglect the periodic structure of this po- 
tential^^ and concentrate our attention only on the three vacua around the origin, 
here denoted by | ) (the one near the origin) and | ±1 ) (the other two). Around 
the vacuum | ), the admitted quantum kink states are 



L) =\ Kqi ) and | S ) 



0,-1 



) 



together with the corresponding antikink states | L) =| Kifi) and | S") =| K^i^). 
The topological charges of these kinks are different, and given by 
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The classical masses of the large and small kink get splitted when we switch on A, 
and their exact value can be easily computed 
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The expansion of this formula for small A is given by 



L,S 



A 71 



where the first term is the classical mass of the unperturbed Sine-Gordon kink 



(6.156) 



6.1 The embarassment of the riches 

Let's now apply the semiclassical formula ()2.12j) for obtaining the spectrum of the 
neutral particles at the vacuum | 0). For the form factor F'^iiO) = /^^(^Tr — 6) of 
the large kink ()6.152|) (with A < 4/i) we have 



A-K 1 sinh [a ^ (ot - 9)] 
~]3 171-9 sinh vr ^ {i7i ~ 9 

' L m ^ 



(6.157) 



^^The periodicity of the potential obviously implies that the following analysis applied as well 
to any other vacuum of this theory. 
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where 



a 



2 arctan 



'4/i + A 



4/i - A ' 

while m and Ml are given by ()6.151|) and ()6.155|) . respectively. For the form factor 
^ss^^^ ~ /|'s(^^ ~ ^) °f small kink ()6.153|) (with A < 4/i) we have instead 

47r 1 sinh [a ^ {ztt - 6)] 



13 m-e sinh [% ^ {iir - 6)] 



(6.158) 



where 



a 



2 arctan 



'4/i - A 



4/i + A ' 

while m and Ms are given by ()6.151|) and ()6.155|) . respectively. 

By looking at the poles of these expressions within the physical strip, it seems 
that there are two towers of neutral particles at the vacuum | ) : the one coming 
from the bound states of the \ LL) 
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2Ml sin ul 
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KTlL < TT- 
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(6.159) 



the other coming from the bound states of | 5 S* ) 
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(6.160) 



As a matter of fact, this situation is not peculiar of the Double Sine Gordon model 
but it occurs each time that there are kinks of different masses originating from 
the same vacuum. Consider, for instance, a simplified version of a three vacua 
configuration, realised by the potential (Figure | 



(^ + 2hf 
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(p > a . 



(6.161) 



The explicit configurations of the long and short kink of this potential are pretty 
simple 
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(6.162) 



and their classical masses are 
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(6.163) 
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Figure 23: Shape of V{ip) with a = b (left hand side) and with a > b (right hand 
side). 



The form factors can be easily computed 
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By looking at the pole in the physical strip of the above amplitudes, it seems then 
that there are two particles on the vacuum | ) , whose masses are expressed by the 
formulas 



niL = 2m a sin 



TT 



ms 



2m b sin 



TT 



(6.166) 



2a2 ' " 262 • 

When a = b, the two masses coincide but, as we discussed in the introduction, the 
corresponding state cannot be degenerate. Hence, one of the two spectra in ()6.166|) 
has to be spurious. The same conclusion applies, as well, to the mass formulas 
(|6.159|) and (|6.16(J|) of the Double Sine Gordon model^'^. But, what went wrong in 
this case with the semiclassical formula? 

To understand the origin of this discrepancy, notice that each kink solution knows 
only half of the shape of the vacuum state from which it originates: for instance, 
the long kink starts its "motion" from the minimum at the origin, but its next 
values are determined only by the shape of the potential on its right. As far as this 
kink solution is concerned, the shape of the potential to the left of the origin could 
be arbitrarly changed without effecting the behavior of this solution. The same 



^•^The argument concerning the non-degeneracy of the neutral states invahdates the conclusions 
previously reached on the spectrum of the Double Sine Gordon model ^5] , and it makes somehow 
obvious the finding of the numerical analysis on this model presented in ^2! ■ 
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considerations also apply to the short kink, which is determined only by the shape 
of the potential on the left of the origin. 

The above observation means that, when we employ the long- kink solution to 
extract the mass spectrum, it is as we are referring to a potential which is not the 
actual one. It is rather a potential Ul{<^), whose values for (yj < are obtained by 
the specular image of those for > of the original potential. Viceversa, when we 
employ the short-kink solution, it is as we are referring to a potential Ufi{ip), whose 
values for (f > are the specular image of those for (f < 0, which determine the 
short-kink solution. For the long and short kinks of the above example, the fictitious 
potentials VL^f) and Vs{f) reconstructed by the semiclassical solutions are the ones 
shown in Figure EU Similar fictitious potentials can be extracted, as well, for the 
Double Sine-Gordon model. Hence, no wonder that employing eq. ()2.12|1 in the case 
of kinks with different mass, each of them gives rise to a different spectrum of bound 
states on the same vacuum. 

Saying the things differently, at the leading order in the coupling constant in 
which the semiclassical form factor ()2.12|1 was computed, the short and long kink 
states are invisible each other. They start to be aware of the existence of the 
other only at the next leading order in the coupling constant. For instance, in the 
expression (j2.25p involving the long-kink, the first subleading terms are given by the 
matrix elements with a couple of short kink and antikink state, as 

{L{ei) I (^(0) \ LSS) {LSS I (^(0) I L)...{L\ (^(0) | 1(62)) . (6.167) 

With Ml > Ms, these terms (as well as all the others, obtained by insering more 
times the couples SS) are always present, no matter which are the values of the 
external rapidities 61 and 62 of the long kink. In this case, it becomes then rather 
artificial to pin down their presence by simply appealing to the perturbation ex- 
pansion in the coupling constant. Their presence, however, spoil the possibility of 
obtaining a close differential equation for the form factors of the kink, as the one 
that has led to the semiclassical expression (j2.12|) . In principle, they can be taken 
in acccount by employing the path integral formalism discussed, for instance, in [23] 
but, in practise, this can be a rather paintful and ferociously complicated procedure. 
So, for the purpose of this paper, much better to content ourselves with the possi- 
bility of identifying the mass spectrum according to the heuristic considerations of 
the next section. 
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(A) 



(B) 



Figure 24: Fictitious potentials V^lip) (A) and Vs{ip) (B) The dashed line gives the 
values of the original potential in the regions not seen by the each kink. 

6.2 The importance of being small 

Once it has been clarified the origin of the discrepancy of the spectra coming from the 
two kink solutions, it remains to understand what is the correct spectrum of bound 
states. Although the exact expression of the mass formula has remained elusive to us, 
we would like to show that the spectrum can be studied, in a relatively simple way, 
at least in two different cases: (a) when the asymmetric kinks have approximately 
the same mass; (b) when the mass of one of them is much smaller than the mass of 
the other. 

Let's consider first the ). This situation can be realised starting by a 

symmetric configuration of the potential (which we assume to be an even function 
V{(p) = V{—(p)) and slightly deforming it by an infinitesimal deformation X6V{(f), 
with 6V{(f) odd under if — ^ —if. Switching on A, the effect of the new term is to 
decrease one the maxima of the potential and to increase the other. This is, for 
instance, what happens in the Double Sine-Gordon for small value of A. Under 
this deformation, the masses of the kinks changes as Ml^s — M ± AA^. Denote 
by I 6)l and by | b)s the bound state of the long and the short kinks which, in 
the unperturbed theory, have equal mass. The actual breather of the unperturbed 
theory is a linear combination of these two "half-breathers" | &)l,5 ^ a combination 
that can be determined by imposing that the state is an eigenvector of the parity 
transformation P (P^ = 1). Let's assume, for instance, that the bound state is 
expressed by the combination 




(6.168) 



In the unperturbed potential, its mass can be equivalently written as 



= - 



{niL + ms) . 



(6.169) 
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Switching now SV, at the first order in A the state ()6.168p does not change. The 

masses and m/j receive, instead, a hnear correction of opposite sign, 

/ /m\m /m\\ 
5l5 ^ ± 2sin — - — cos — . 6.170 
V \2MJ M \2MJJ ^ ' 

Plugging these corrections into ()6.169p . the mass of the breather remains then un- 
changed. This result matches with the first order Form Factor Perturbation Theory 
|161 117j . as it can be seen by employing the parity operator P 

\{B\5V{^) \B) = X{B\P{P6V{^)P) P\B) = -\{B\5V{^) |5) = . 

(6.171) 

Let's consider now the case (b), i.e. when one of the kink is much heavier than 
the other. We would like to present a series of arguments in favour of the thesis 
that, in this circumstance, the semiclassical spectrum is essentially determined by 
the short-kink solution, i.e. 

^ 2M^sin [us ^) , < < vr^ , (6.172) 

with the stable part of the spectrum obtained only for ng < 2. The arguments are 
the following 

1. The actual mass of the bound state depends on both Ms and Ml. However, 
repeating the argument presented in Section IH.H the correction induced by 
the long kink is expected to be suppressed with respect to the one of the short 
kink by the ratio j • Therefore, making heavier the mass of the long kink, 
it can be forseen that its influence on the mass of the neutral particle should 
becomes less and less relevant. In particular, when Ml oo, the long kink 
decouples from the theory and the masses of the breathers are only determined 
by the dynamics of the remaining short kink. 

2. Another argument in favour of the formula ()6.172j) directly comes from the 
Watson equations satisfied by the form factors of the short and long kink 
solutions. Let's first consider the scattering process of the small kink-antikink 
state I S{9i)S{92) )■ As far as we have the inequality 

Ms cosh^ < Ml cosh^ , (6.173) 

in the center of mass (defined by | = 9i = —62 and | = 6'i = —O^)^ there is no 
possibility of converting the above state into a | L{9i)L{92) ). Therefore the 
scattering process of this short kink-antikink state can only be elastic 

I 5(^1) 5(4)) I 5(4)5(4)) . (6.174) 
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Figure 25: Elastic scattering amplitude of the short kink-antikink state. 

The elastic range of 6 obviously enlarges by making the mass of the long kink 
heavier. In this region, the Watson equation which is satisfied by any form 
factor of the short kink-antikink state becomes then 

Fg{e) = Sssie)Fgi-e) , (6.175) 

where 6*55(6') is the elastic S*- matrix of the process ()6.174p . Assuming that 
this elastic process leads through the exchange of the scalar particles, we have 
the diagram of Figure ESI 

Eg. (jfi.l75|l implies that the ratio F^g{9) / Fgg{—9) is a pure phase for the 
real values of 6 below the inelastic threshold given by eq. ()fi.l73j) . perfectly in 
agreement with the semi classical result. If we now trust the semi classical result 
of the form factors also for complex values of the rapidity (in particular for 
those concerning the location of their poles), we arrive to the mass spectrum 

mm- 

Consider now the long kink-antikink scattering. First of all, notice that 
for the non-integrability of the theory, the scattering processes of the state 
I L{9i) L{62) ) always involve, as a final state, | 5'(^i) S{92), ) 

\m)L{92)) ^ \S{9^)Si92)) . (6.176) 

i.e. the production process is always present. 
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Figure 26: Elastic and production amplitudes of the long kink-antikink scatterings. 

Going to the center of mass, the rapidity of the outcoming small kink-antikink 
state is determined by 

Ms cosh^ = Ml cosh^ . (6.177) 

In addition to the production process, in the scattering of the long kink- 
antikink state, there is also its elastic part 

\L{e,)L{e2)) — |L(^i)Z(^2)) . (6.178) 

For values of 6 below other inelastic thresholds, the Watson equation satisfied 
by the form factors of the long kink-antikink state is then 

Ffz(^) = SLLie) F^i^i-e) + sue) Fg{-e) , (6.179) 

where Sll{G) is the S'-matrix relative to the process ()6.178|1 . whereas Slr{0) 
is the one of ()6.176j) . Assuming that, also in this case, the scattering pro- 
cesses are dominated by the exchange of the scalar particles, we have the 
diagrams of Figure EHl In contrast with the short kink case, this time the 
ratio F^j^{9) / F^^{—9) can never be a pure phase, not even for real values of 
6. Hence the semi-classical result cannot be the correct one: in fact, from the 
ratio of the long kink form factors, one always gets a pure phase expression. 

Notice that in the production diagram | LL) ^\ SS) it appears the same 
3-particle coupling than in the elastic scattering | SS) ^\ SS). This implies 
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Figure 27: Relation between the resonance angles of the short and long kink-antikink 
state. 

that the intermediate particles of the two processes are the same. Assuming 
that their masses are those extracted by the short kink-antikink bound states, 
we can predict where the correct position of the poles of the long kink-antikink 
form factors should be: denoting the position of these poles by iu'g^ = Z7r(l — 
C§s)^ with = n-^, and = Z7r(l — ■C^z)' ^^e resonance value 
determined by the relation 



shown in Figure ()27p. It is easy to see that this relation is similar to the law 
of refraction of hght between two media of refraction indices and M/j. 

3. In order to extract the spectrum of the neutral particles at a given vacuum \ a), 
the actual thing to do is, of course, to quantize the time- dependent solution 
nearby the corresponding minimum ip'^a^ of the potential. To this aim, one 
initially needs to solve the equation of motion 




(6.180) 



dU 




dt"^ dx"^ 



dip ' 



by requiring the following properties of the solution: 



• To be a periodic function in time, with a certain frequency uj. 

• To be localised around the minimum (p^a^ of the potential. 
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To have a finite value of its energy, given by 



E 



1 f d(p\ i f dip 



2\dt 2\dx 



dx . (6.182) 



Tlie finiteness of the energy implies that, at x ^ ±00, the solution should 
necessarily go to ip^a \ where the potential vanishes. However, we must 
also require that E < 2Ms, which is the lowest energy threshold in the 
neutral sector. Without this condition, in fact, the asymptotical time 
evolution of the solution will consist of a state of small kink and small 
antikink, moving a part with respect each other, with a velocity fixed by 
the excess of energy with respect to their rest mass. This is the classical 
equivalence of a decay process. 

Having stated these conditions, the problem of solving eq. ()6.181|) is equivalent 
to find the small oscillation of a string in the funnel of the potential U{ip) 
relative to the minimum ipii \ with its ends (at x = ±00) kept fixed at the 
value of its bottom fFigure I^Hj). During its swinging, the string explores and 
probes the shape of the potential nearby the valley of the minimum y^i^^ The 
oscillations, however, cannot be too wide, otherwise the motion will violate 
the energy condition E < 2Ms- Notice that, in this analogy, the static kink 
solutions are nothing else than strings which, at x = — 00 are at the bottom 
of one valley, whereas at x = cxo, after overpassing one of the peaks of the 
potential landscape, are in a neighbouring one. 



In order to solve ()6.18H) . let's pose (p{x, t) = (pa + ^{x, t) and expand corre- 
spondingly the potential near this minimum 

^(^r+r/) = IujW + yv' + ^v' + ^v' + --- . (6.183) 

Notice that the asymmetry of the potential with respect to its two barriers 
(i.e. the mere fact that there exist the small and the long kink) is encoded 
in the non-zero coefficients of the odd powers of 77. Any smooth, real-valued 
solution of ()6.181|) that is periodic in time with frequency u > has a Fourier 
representation 

r]{x,t) = ^^a„(x) exp{inut) , a_„(x) = a*(x) . (6.184) 

All the coefficients a„(x) must be localised, i.e. they should vanish when 
X ±00. Substituting this expansion into ()6.181|1 and using eq. ()6.183j) . one 
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Figure 28: Landscape of the potential seen by a string oscillating in the central valley, 
relative to the minimum ip^a^ . 



reaches the non-hnear equations for a„'s 



/ o2 , 2 2 2\ 
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k m 



(6.185) 



Spatially uniform infinitesimal solution of ()6.181|) oscillate with frequency uq, 
and it is possible to prove that ()6.18H) admits no breather solution for u > uq 
|26j . Classically, this result is expressing the fact that the actual mass is always 
smaller than the one fixed by the curvature of the potential. Hence, one can 
focus on the interval < u; < ujq. In order to find a small parameter e to make 
a reasonable expansion, previous studies of this equation^^ (see, for instance, 
0123 EH!) suggest to use 



and to introduce a slightly modified asymptotic expansion in terms of the 

^^It is interesting to observe that in |27ll28j . the authors actually argue about the •non-existence, 
stricly speaking, of a stable solution of eq. Ht).181|) . However, the stability breaking of the solution 
occurs through exponentially small terms, invisible to any perturbative orders, which correspond 
to a radiation rate of the perturbative solution. Due to the extremely slow decay of these processes, 
for any practical purpose, one can safely ignore this mathematical subtlety, as also admitted by 
the same authors. 




(6.186) 
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dimensionless variables: 

The solution of (j6.185p will be expressed in terms of a series expansion in e of 
all the terms a„. For instance, in the case of Sine-Gordon model, the breather 
mode (at the leading order in e) is given by [iHl EHj 



4 / esinfco-ot/vT+e^) , 
m(x^t) = -arctan ^ " ^ \ ^ 6.188 

^ \cosh(ec^oa^/v/rT^) 

4 

~ — sech,^ sin[r(l — e^/2)] + ■ ■ ■ . 

Once such solution of ()fi.l85p has been found, the next steps will be to compute 
its classical energy E and its action 

along one period. Referring once again to the Sine-Gordon model, the results 
are 

3 



E = 2M(e- — + ■■■) , (6.190) 



47iM f e3 

W = + 

m \ 6 

where M is the mass of the soliton, alias the heights (per unit length) of the 
barriers of the potential landspace. Eliminating now e, one has 

E^,„(}^-L(I^y-(K)\...] , (6.191) 

\^27r 24 VM/ \27i J J ^ ' 

The familiar mass spectrum of the Sine-Gordon model is finally obtained by 
imposing the quantization condition 

W = 27rn . (6.192) 

If the same calculations are repeated for ip^ theory |2H] , one finds 

m 



—= sech,^ cos r — sech^,^ 



-|- — sech i cos 2r H -= sech ^ cos 3r -|- 

3 6v^ 



(6.193) 
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with the classical energy and the action of the swinging string in the valley of 
this potential given by 



As in the Sine-Gordon model, also in this case the energy is expressed in 
terms of the mass of the soliton, alias in terms of the height of the barrier of 
the potential. Eliminating e and imposing the quantizion condition ()6.192p . 
one arrives, as before, to the series expansion of the usual formula = 



The two examples above should have made clear that, in the general case, the 
energy and the action of the classical solution will be expressed in terms of the 
lowest heights of the potential, a feature which is pretty intuitive. In the case 
in which these heights are the same (like the Sine-Gordon case), there is only 
one energy scale, given by the mass of the soliton. But, also in the case of y?^, 
there is only one energy scale, given by the height of the potential between 
the two vacua: on the left (or on the right) of each of them, there is in fact an 
infinite barrier. 

The same one-scale situation is expected to be valid when the potential has a 
barrier much higher than the other: in this case, the small oscillations of the 
string will be essentially determined by the barrier given by the lowest peak 
(the other barrier wU only induce corrections at higher order in e). Correspond- 
ingly, the spectrum of the neutral bound states should essentially coincide with 
the one extracted by the lowest kink. 

This discussion should also clarify the reason of the difficulty to find an exact 
expression of the mass spectrum when the kinks have a comparable mass. In 
this case, in fact, there are two scales in the problem, i.e. Ml and Ms, and 
the energy together with the action of the string will be non-trivial functions 
thereof. To find their expression, at least in a concrete example, is an inter- 
esting open problem on which we hope to come back in the future. From a 
practical point of view, notice that, when Ml is not much larger than Ms, the 
masses obtained by employing either eq. ()6.159l or eq. ()6.160|) are always very 
close each other and they can provide an indication on the actual value of the 
masses of the bound states. 

Even though it seems rather difficult, in general, to determine where the 
crossover from the two-scale scenario to the one-scale scenario takes place. 



W 
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(6.194) 



2Msin (n^). 
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when Ml/ Ms > 2 one should be relatively safe by taking Ms as the only scale 
of the problem^^. 

7 The BLLG potentials 

In this section we will briefly discuss the particle content of an interesting class 
of potentials, introduced by Bazeia et al. in [21]. These potentials, which can be 
expressed in terms of the Chebyshev polynomials of second kind, are closely related 
to the if'^ potential (in its broken phase). In fact, they are obtained from this theory 
by using the so-called deformation procedure, explained in [20] . In the following we 
denote them as BLLG potentials. In terms of the dimensionless coordinates and the 
dimensionless field previously used (x^ — > mx^, cp "^V^)? the Lagrangian of the 
BLLG models is given by 



4 

c = - 



1,„ ,2 



2 



(7.195) 



where 



Ua{^) = ^(l-^')'K-iM , (7.196) 



sin[(a + 1) arccosv?] 
sin(arccos (/?) 



The potentials of these theories fall in two classes, depending whether a is an odd 
or an even integer: those with a odd are like (p"^ in its broken phase, i.e. with a 
maximum at the origin, whereas those with a even are like (f^ potential, with a 
zero at the origin. Some of these potentials are drawn in Figure EHl Following the 
analysis done in [211, ^^^'^ briefly summarise their main properties. 
The minima of the above potentials are localised at 

(^(°) = cos ( ^— ^vr ) , A;= l,2,...,2a + l . (7.197) 



By taking into account the periodicities ipf^ = — </?^°^2-fc '^a'+k ~ ^'a+2-ki 
easy to see that the number of distinct (positive) values are a + 1, all negative values 
obtained by reflection. For a even, one has this set of zeros 



^^The string oscillations are localised in a region order e. It is hard to immagine a string 
oscillating twice higher than the first barrier, whose height for unit length is expressed by Ms, still 
keeping its energy lower than 2Mg. 
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(A) (B) 
Figure 29: BLLG potentials Usiff) (A) and Ug{(^) (B). 

whereas for a odd 

= 1, , . . . , ^fJX , -^it , . . . , -^f, -^T = -1} . (7.199) 

These minima are the vacuum states | A; ) (A; = 1, 2, . . . , a + 1) of the corresponding 
quantum theory. 

The nice thing about these potentials is that, thanks to the deformation proce- 
dure, all kink solutions are explicitly known. They are expressed as 

^Jx) = cos-^^—^ — , 7.200 

a 

where k = 1,2, ... ,a and ip{x) G [0, vr] is the principal determination of arccos tanh(a:). 
These topological configurations interpolate between the vacua ^fli (reached at 
X = —oo) and ip^^^ (reached a.t x = oo). The solutions ipk+a{x) are the correspond- 
ing anti-kink configurations. In order to compute their mass, it is useful to introduce 
the so-called super-potential Wa{(p) and write the potential Ua{(p) as 

1 fdW„^ ' 



U.i,)--,[-^) . (7.201) 

This is always possibile since Ua{(p) > 0. For a ^ 2, the explicit expression of Wa{<p) 
is 

Wai^p) = ] [(a'(l - <^') - 2) cos(a arccos (^) (7.202) 
a^(a^ — A) 

—2aip \/l — Lf'^ sin (a arccos yj) ] , 

whereas, for a = 2, 

W2{^) = ^<^'(2-(^2) . (7.203) 
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(A) (B) 
Figure 30: Masses of the kinks in the BLLG potentials U^{ip) (A) and Ug{ip) (B). 



In terms of the super-potential, the mass of the /c-th kink is given by 

i3 



16 



T 



k+l) 



(7.204) 



By using (j7.203p . we arrive to the final expression of the kink masses 



A a2(a2-4) 



sm 



{k - l)-K 



sm 



a ) 



, a^2 

(7.205) 

whereas Mi = ^m^/X for a = 2. 

The kink with the lower mass is the one connecting the farthest vacua of the 
potential, either on the left or on the right. The kinks with the higher mass are, 
instead, the ones related to the vacua nearby the origin. Since the masses of the 
kinks are proportional to the heights of the potential, by making an istogram of the 
their values - each of them placed at the middle of the two vacua interpolated by the 
corresponding kink - one expects to get a scheleton form of the original potential, 
as it is indeed the case, see Figure EDI 

Let's discuss now the spectrum of the neutral particles at each vacuum | k). 
Notice that all vacua, except the most external ones, have the same curvature. More 
precisely, reintroducing the correct dimensional units and denoting the curvature by 
u"^, one has 



Co", 



a+l 



4m 



UJi. 



m 



(7.206) 



The perturbative mass of the neutral particles at the vacua | 1 ) and | a + l) is 
twice larger than the one at the other vacua. Together with the lowest value of the 
mass of the kink that originates from the most external minima, it is clear than the 
most sensitive situation for the existence of the neutral particles happens at the two 
farthest vacua. 



16 



In this formula we have restored the original dimensional quantities. 
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Referring to the discussion of the previous section, as far as the mass of the large 
kink is much higher than the mass of the short kink, the spectrum of the neutral 
particles at the corresponding vacuum is determined by the smallest one. If these 
masses are comparable, if one prefers, can use instead their average. The conditions 
on the relative weigths of the masses can be easily checked for the BLLG potential. 
Since, in both cases, it is always the lowest mass that matters, we can simplify the 
notation, at least, by always employing the masses of the smaller kinks. To have 
neutral particles at the vacuum | /c ), one has then to check the condition 

6 = ^ < vr , (7.207) 

where is the mass of the small kink (or antikink) which has the vacuum | k ) 
as asymptotic limit at x —oo. Specializing the above formula to the vacua | 1 ) 
and I a + 1 ) , one has 

2a^(a^-4) ^ 
4 - sin^ - 



6 = ^^^^ ■ ("08) 



Concerning the values of at the other vacua, notice that for all vacua \k) io the 
left of the origin, the small kink are always the antikink of </3fc(x). For those to the 
right of the origin, the small kink is instead the kink ipk{x) itself. We have 



= , ""Yu TT^A ' k = 2,...ra (7.209) 

4_a2fsin2(^ + sin2(^'l 



a 

a'(«'-4) A , ~ , 

Ek = 7 ^ ;7 , k = a + 1, 



4 - a"^ sm — + sm 

\ a a 

where a = | + 1 for a even, whereas a = for a odd. For instance, in the case of 
the potential Ug^ip), with the notation ^ = E,m? /{jiX), one finds 

= 725.2.. 

6 = 362.6.. 

^3 = 50.98.. 

^4 = 22.06.. 

(7.210) 

^6 = 14.67.. 

= 22.06.. 

is = 50.98.. 

= 362.6.. 

U = 725.2.. 
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Hence, there are a series of nested equations relative to the bound states on the 
various vacua: for instance, if 

A < i^' ' 

there is one particle on the vacua | 1 ) and | 10 ) and two particles on all the others. 
Increasing the value of ^ to the interval 

fr^ < A < , 

the farthest vacua do not have neutral excitations any longer, the vacua | 2 ) and 
I 9 ) have one particle, while all the other have two bound states. Increasing ^, there 
is a progressive emptying of particles on the various vacua, so that, when 

- > 



there are no more neutral particles on all vacua. 



8 Conclusions 

In this paper we have used simple arguments of the semi-classical analysis to in- 
vestigate the spectrum of neutral particles in a quantum field theory with kink 
excitations. Leaving apart the exact values of the quantities extracted by the semi- 
classical methods, it is perhaps more important to underline some general features 
which have emerged through this analysis. One of them concerns, for instance, the 
existence of a critical value of the coupling constant, beyond which there are no 
neutral bound states. Another result is about the maximum number n < 2 of neu- 
tral particles living on a generica vacuum of a non-integrable theory. An additional 
aspect is the role played by the asymmetric vacua and by the asymmetric kinks. 

There are several interesting open problems which deserve a further investigation. 
An important open question is to find the exact mass formula (if it exists) when the 
asymmetric kinks have a comparable value of their masses. This goes together with 
the problem of finding a convenient way of taking into account higher order terms 
in the form factor expression of the kinks. Another challenging aspect concerns 
the refinement of the analysis of the resonances and of the corresponding decay 
processes. Under this respect, it may be possible that useful insights will come in 
the future from a numerical analysis. 
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